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ABSTRACT 


Motivated  by  the  sizable  ledger  of  ships  sent  to  the  ocean  floor  without  ever 
sustaining  a  direct  hit  during  World  War  II,  a  heightened  interest  in  ship  shock 
survivability  spread  throughout  the  Naval  Engineering  community.  As  a  result,  over  the 
last  fifty  years,  Live  Fire  Test  &  Evaluations,  otherwise  known  as  ship  shock  trials,  have 
been  conducted  in  order  to  determine  the  seaworthiness  of  each  new  class  of  ship 
commissioned  in  the  U.S.  Fleet.  While  beneficial  in  determining  the  overall  survivability 
of  a  ship  and  its  mission  essential  equipment  in  a  severe  shock  environment,  these  Navy- 
mandated  tests  pose  serious  danger  to  the  crew,  ship  and  environment. 

As  an  alternative  to  these  labor  intensive,  costly  and  time  consuming  at-sea  tests, 
the  recent  advances  in  computer  processing  power  have  made  it  possible  to  employ  finite 
element  methods  involving  complex  geometries  in  the  modeling  and  simulation  of  shock 
response  for  the  ship  and  surrounding  fluid.  This  thesis  examines  the  accuracy  of  shock 
simulation  predictions  as  compared  to  the  ship  shock  trials  conducted  on  USS  WINSTON 
S.  CHURCHILL  (DDG-81).  An  investigation  of  the  effects  of  sensor  location,  damping 
and  shot  geometry  is  presented  as  validation  of  the  Naval  Postgraduate  School  modeling 
and  simulation  methodology. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

To  the  dismay  of  many  World  War  II  naval  leaders,  scores  of  ships,  both 
combatants  and  those  in  service  with  the  merchant  marine,  were  sent  to  the  ocean  floor 
without  ever  taking  a  direct  hit  from  an  enemy  bomb,  mine  or  torpedo.  These  ships  sank 
as  a  result  of  underwater  explosions  (UNDEX)  occurring  in  the  surrounding  waters 
adjacent  to  the  ship.  In  the  ensuing  years  a  buzz-phrase  sharing  the  title  of  a  training  film 
on  underwater  explosions,  “near  is  close  enough”,  proved  to  signal  a  shift  in  naval 
warfare.  With  the  ability  to  effectively  deliver  ever-increasing  charge  sizes  it  became 
readily  evident  that  hitting  the  hull  of  the  ship  was  no  longer  as  important  as  once  had 
been  the  case.  It  had  been  generally  accepted  throughout  the  Fleet  that  one  of  the  best 
ways  to  sink  a  ship  was  to  open  a  large  hole  in  the  hull  beneath  the  waterline.  This  was 
done  with  the  expectation  that  the  ensuing  flooding  would  reduce  stability  to  a  point  from 
which  the  ship  could  not  recover.  Delivering  a  direct  hit  to  a  weapons  magazine  or  fuel 
storage  tanks  that  would  facilitate  the  cascade  of  internal  explosions  and  ultimate 
catastrophic  loss  due  to  conflagration  was  also  deemed  highly  desirable  and  practical. 
However,  through  insightful  analysis  of  the  wartime  ship  losses  suffered  during  the  first 
half  of  the  20th  century,  it  was  determine  that  incident  shock  wave  and  bubble  pulse 
forces  resulting  from  UNDEX  events  were  one  of  the  primary  initiators  of  structural 
damage,  material  failure  and  ultimate  loss  mechanism  in  the  sinking  of  numerous  ships 
[Ref.  1], 

Over  the  last  fifty  years  much  research  has  been  accomplished  in  the  UNDEX 
field,  resulting  in  a  greater  appreciation  of  the  true  power  encapsulated  in  the  UNDEX 
shock  phenomena.  Accordingly,  having  understood  the  necessity  for  ships  that  were 
resilient  in  an  UNDEX  environment,  specifications  were  established  for  the  design  and 
testing  requirements  of  all  naval  surface  combatants.  The  Department  of  the  Navy  set 
forth  guidance  for  shock  hardening  of  surface  ships  in  OPNAVINST  9072.2  [Ref.  2], 
with  additional  requirements  delineated  in  NAVSEA  0908-LP-000-3010A  [Ref.  3]  and 
MIL-S-901D  [Ref.  4],  Carried  out  in  the  summer  of  2001,  the  DDG-81  Ship  Shock 
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Trials  are  the  latest  set  of  Live  Fire  Testing  &  Evaluations  (LFT&E)  to  be  conducted  in 
fulfillment  of  these  standing  requirements. 

Referred  to  as  “shots”,  a  series  of  underwater  explosions,  produced  by  the 
detonation  of  charges  placed  at  varying  locations  in  the  water  adjacent  the  ship,  are 
designed  to  test  the  ship  at  “near  combat  conditions”  [Ref.  2],  The  response  of  the  ship, 
weapons  systems,  specific  equipment  and  the  crew  are  all  measured  and  recorded  in  order 
to  evaluate  their  performance  in  a  shock  environment.  Taken  from  one  of  the  first  ships 
in  the  class,  or  from  a  ship  incorporating  major  design  changes  during  construction,  this 
data  is  then  analyzed  and  recommendations  are  made  for  the  alteration  of  existing  ships 
or  for  a  change  in  the  design  of  subsequent  ships  to  be  built  within  that  same  ship  class. 

Even  though  these  ship  shock  trials  provide  a  true  record  of  the  system  response 
of  the  ship  as  well  as  excellent  training  for  the  ship’s  crew,  they  are  very  costly,  and 
inherently  dangerous.  Such  events  require  extensive  planning  and  coordination  and  are 
potentially  damaging  to  the  ship  structure,  electronics  and  multi-million  dollar  weapons 
systems.  Inasmuch  as  these  shock  trials  are  good  measures  of  the  ship’s  potential 
performance  in  a  shock  environment,  they  are  limited  by  the  safety  risk  involved  and  thus 
only  test  to  two-thirds  the  design  limit.  These  limitations  of  the  LFT&E  program  raise 
concerns  over  the  validity  of  the  ship  shock  trials  and  their  associated  costs,  which  could 
range  as  high  as  5%  of  the  $950million  delivery  cost,  as  in  the  case  of  the  USS  JOHN 
PAUL  JONES  (DDG-53)  ship  shock  trials  conducted  in  1994  [Ref.  5]. 

In  recent  years,  unprecedented  advances  in  computer  modeling  and  simulation 
have  created  the  potential  to  mitigate  some  of  the  costs  associated  with  the  LFT&E 
activities  through  the  use  of  virtual  shock  environment  analysis  [Ref.  6].  By 
implementing  these  current  technologies,  simulations  that  accurately  predict  the  initial 
peak  response  of  a  surface  ship  subjected  to  an  underwater  shock  event  enhance 
traditional  analysis  methods  and  hold  great  promise  in  replacing  certain  types  of  at  sea 
live  fire  testing  [Ref.  7].  With  the  two  major  elements  of  the  modeling  and  simulation 
process,  model  refinement  and  computer  runtime,  amounting  to  only  a  fraction  of  the 
ship  shock  trial  costs,  engineers  can  subject  the  finite  element  ship  model  and 
corresponding  fluid  mesh  to  an  exhaustive  battery  of  simulations  over  an  extensive  range 


2 


of  charge  sizes  and  geometries.  These  simulations  conducted  in  the  virtual  UNDEX 
environment  allow  for  evaluation  of  the  ship  system  at  and  beyond  its  design  limits, 
without  bringing  harm  to  the  ship,  crew  or  the  environment.  This  approach  has  the  added 
benefit  of  incorporating  predictive  results  obtained  from  the  simulations  into  the  final 
stages  of  the  ship  design  spiral.  Making  corrective  changes  while  the  ship  is  still  in  the 
construction  phase,  rather  than  during  the  post-production  timeframe  as  is  done  with 
current  LFT&E  results,  greatly  reduces  the  rework  costs. 

Though  not  considered  reliable  enough  at  this  time  to  completely  replace  the 
LFT&E  process,  it  is  expected  that  the  use  of  simulated  UNDEX  events  will  continue  to 
be  used  as  a  predictive  design  tool.  The  insight  gained  in  the  virtual  UNDEX 
environment  would  eliminate  the  need  for  broad  scope  shots  and  foster  focused 
investigation  of  UNDEX  events  through  use  of  scalable  charges  placed  at  specific 
locations  corresponding  to  the  points  of  interest  found  in  pre-shock  trial  simulations. 


B.  SCOPE  OF  RESEARCH 

Using  the  data  obtained  from  the  shock  trials  conducted  on  USS  WINSTON  S. 
CHURCHILL  (DDG-81)  in  June  of  2001  as  a  basis,  this  paper  serves  as  further 
validation  of  the  modeling  and  simulation  methodology  established  by  the  Shock  and 
Vibration  Computation  Laboratory  at  the  Naval  Postgraduate  School  (NPS).  Recent 
work  completed  by  Schneider  [Ref.  8]  at  NPS  concentrated  mainly  on  the  refinement  of 
the  modeling  and  simulation  process  of  the  DDG-81  as  well  as  analysis  on  the  overall 
vertical  response  velocities  at  the  shipwide  level.  This  paper  will  expound  upon  the  work 
presented  by  Schneider  in  June  2003.  Using  the  NPS  modeling  and  simulation  process, 
this  paper  further  investigates  the  results  of  localized  vertical  velocity  response  in  the 
critical  Combat  Information  Center  (CIC)  area  as  well  as  the  shipwide  velocity  response 
in  the  athwartship  direction.  Comparisons  between  the  measured  ship  shock  trial  data  and 
the  simulated  response  predictions  were  also  conducted  for  two  sets  of  ship  system 
proportional  damping  coefficients.  The  findings  of  these  comparisons  will  also  be 
presented  herein. 
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II.  UNDERWATER  EXPLOSION  PHENOMENA 


In  order  to  truly  understand  the  devastating  effects  that  are  associated  with  the 
underwater  shock  phenomena  it  is  necessary  to  start  with  some  background  information 
on  this  subject.  Since  there  are  many  complex  layers  to  the  underwater  shock  phenomena 
and  its  corresponding  system  response,  only  the  most  relevant  factors  will  be  presented  in 
order  to  provide  some  degree  of  familiarization. 


A.  FLUID  BEHAVIOR 

There  is  a  defined  sequence  of  events  that  makes  up  the  underwater  explosion 
(UNDEX).  We  shall  begin  with  a  simple  discussion  of  the  underwater  shock  and  where 
and  how  it  comes  exist.  To  set  the  process  in  motion,  a  sudden  chemical  reaction  within 
the  initiator  charge,  typically  composed  of  a  primer  explosive  such  as  mercury  fulminate, 
results  from  an  electronic  or  mechanical  detonation.  This  action  then  ignites  the  high 
explosive.  Commonly  used  high  explosives  are  HBX-1,  RDX,  TNT  and  PETN.  Their 
detonation  and  ensuing  conflagration  causes  a  high  temperature,  high  pressure  gas  to  be 
formed.  Almost  immediately  a  shock  wave  propagates  outward  from  the  nucleus  of  the 
charge  at  a  velocity  on  the  order  of  25,000  ft/sec  [Ref.  9].  With  the  reaction  initiated,  a 
pressure  wave  proceeds  to  moves  through  the  surrounding  explosive  material,  creating 
additional  pressure  waves.  The  generated  explosive  energy  exists  in  a  gaseous  state  with 
temperatures  and  pressures  approaching  3000  degrees  Celsius  and  50,000  atmospheres, 
respectively  [Ref.  10].  This  initial  process  takes  only  nanoseconds  to  occur  in  most  high 
explosives  [Ref.  11].  With  the  pressure  wave  velocity  exceeding  the  acoustic  velocity  of 
the  explosive  material  by  anywhere  from  three  to  fives  times,  a  shock  wave  is  formed. 
This  combination  of  extremely  high  heat  and  compressive  pressures  facilitates  the  self- 
perpetuating  nature  of  the  explosive  process.  The  resulting  shock  wave  is  then  released 
into  the  surrounding  fluid. 

Typically  when  dealing  with  water  in  engineering  applications,  it  is  taken  to  be  an 
incompressible  fluid  in  all  but  the  most  rare  of  cases.  However,  in  UNDEX  applications, 
the  water  immediately  surrounding  the  explosive  charge  actually  compresses  slightly  as  a 
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result  of  the  extreme  pressure  of  the  shock  wave.  This  compression  results  in  a  high 
pressure  shock  wave  in  the  water  that  in  turn  spreads  outward  from  the  charge  location. 
Though  this  shock  wave  initially  moves  through  the  water  at  velocities  much  greater  than 
the  speed  of  sound,  it  quickly  retards  to  match  the  acoustic  velocity  of  the  water. 
Generally  approximated  as  5000  ft/sec,  the  actual  speed  of  sound  through  water  is 
affected  by  such  factors  as  temperature,  hydrostatic  pressure,  and  salinity  of  the  medium 
[Ref.  10].  For  the  simulations  discussed  herein  a  value  of  5057  ft/sec  is  used  in  all  cases. 

The  pressure  wave  generated  during  the  detonation  process  has  an  incredible 
amount  of  force  driving  it  outward  from  the  charge  center.  For  example,  in  the  case  of 
TNT,  it  is  on  the  order  of  2xl06  lb/in2.  Figure  1  shows  an  example  of  the  pressure  profile 
for  a  TNT  charge  [Ref.  9],  The  initial  shock  wave  shows  a  discontinuous  pattern  of 
exponential  decay  as  the  radial  distance  from  the  detonation  point  increases. 


Figure  1.  Shock  Wave  Profiles  for  300  lb  TNT  Charge  [from  Ref.  9] 


A  series  of  empirical  equations  have  been  formulated  to  characterize  the  shock 
wave  pressure  profile,  P(t)  [Ref.  10].  These  formulae,  Equations  (1)  through  (5),  are 
good  from  10-100  charge  radii  and  for  up  to  one  time  decay  constant  after  the  initial 
detonation. 
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where  the  preceding  variables  are  defined  as: 
ti  =  arrival  time  of  shock  wave  (msec) 
t  =  time  of  interest  (msec) 

Pmax  =  the  peak  magnitude  of  the  pressure  of  the  shock  front  (psi) 
0  =  shock  wave  decay  constant  (msec) 

R  =  standoff  distance,  radial  (ft) 

W  =  weight  of  the  explosive  (lb) 

D  =  charge  depth  (ft) 

Ki,  K2,  K5,  Kg,  Ai,  A2  =  constants  specific  to  explosive  type 
Amax  =  maximum  bubble  radius  (ft) 
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B.  GAS  BUBBLE  OSCILLATION 

After  the  shock  wave  is  produced  it  quickly  expands  radially,  leaving  behind  the 
highly  compressed,  superheated  gases  that  formed  it.  This  spherical  gaseous  bubble 
continues  to  expand  to  relieve  its  pressure  until  the  internal  pressure  falls  below  the 
surrounding  hydrostatic  pressure  of  the  water.  In  the  interim  the  bubble  actually  expands 
past  its  equilibrium  due  to  the  momentum  of  the  expansion,  growing  to  nearly  twice  its 
equilibrium  diameter.  Equation  (5)  is  used  to  calculate  the  maximum  bubble  radius 
[Ref.  9], 

This  equation  illustrates  that  the  maximum  bubble  radius  has  a  one-third  power 
relationship  with  the  charge  weight  and  an  inverse  proportionality  to  the  one-third  power 
of  the  charge  depth.  At  the  instant  that  the  gas  bubble  reaches  its  maximum  diameter, 
there  is  a  sizable  positive  pressure  gradient  between  the  bubble  and  the  encompassing 
water,  which  causes  the  gas  bubble  to  implode  upon  itself.  The  bubble  then  shrinks  down 
to  a  point  where  the  pressure  within  the  bubble  is  high  enough  to  prevent  further  collapse. 
A  negative  pressure  gradient  now  exists  between  the  bubble  and  the  water  that  surrounds 
it.  Once  again  the  bubble  attempts  to  expand  to  its  equilibrium  state  and  reaches  a 
maximum  diameter  smaller  than  the  initial  gas  bubble  diameter,  yet  still  larger  than  its 
expected  equilibrium  point,  though  the  overshoot  is  less  than  in  the  first  case.  Henceforth 
this  oscillatory  process  repeats  until  the  energy  contained  within  the  bubble  is  insufficient 
to  continue  the  cycle  or  the  bubble  has  come  in  close  proximity  with  the  free  surface  of 
the  water,  allowing  the  exhaust  gases  to  freely  vent  to  the  air  above.  The  cyclic 
expansion  and  contraction  of  the  bubble  along  with  its  migration  path  to  the  free  surface 
are  shown  in  Figure  2. 

The  vertical  migration  velocity  experience  by  the  bubble  in  its  ascent  to  the  air- 
water  interface  is  calculated  by  using  Equation(6).  The  vertical  velocity  (U)  is  a  function 
of  the  gas  bubble  radius, 

U=4Ma3(t)dt  (fl/sec)  (6) 

a3(t)J 
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where: 

g  =  gravitational  acceleration  constant 
a  =  gas  bubble  radius 

Even  though  the  gas  bubble  pulse  is  highly  dependent  on  charge  geometry, 
specifically  charge  size  and  detonation  depth,  it  is  important  to  the  simulation  of  the 
UNDEX  event.  This  is  especially  true  in  the  case  of  the  DDG-81  shock  trial  simulations 
since  the  gas  bubble  pulse  has  a  low  oscillation  frequency  that  approaches  the  values  of 
the  first  bending  mode  of  the  ship.  It  could  potentially  result  in  even  more  destructive 
forces  than  the  incident  pressure  wave,  given  the  proximity  of  the  phenomena  [Ref.  10]. 


Shock  Wave  Bubble  Pulse 


Figure  2.  Gas  Bubble  Oscillation  and  Migration  Path  [after  Ref.  9] 
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C.  CAVITATION 

In  general,  cavitation  is  described  as  the  phenomenon  that  occurs  when  there  is  a 
region  of  absolute  negative  pressure  present  in  a  body  of  water  at  a  nearly  constant 
temperature.  Vapor  cavities  are  formed  where  the  negative  pressure  causes  the  tensile 
force  of  the  water,  which  cannot  be  sustained,  to  implode,  in  turn  causing  high  pressure 
pulses  to  emanate.  During  an  UNDEX  event  there  are  two  types  of  cavitation  present, 
bulk  cavitation  and  local  cavitation.  The  first  type,  bulk  cavitation,  can  be  thought  of  as  a 
large  area  of  low  pressure  at  the  free  surface  of  the  water  forming  above  the  charge 
detonation  point,  while  the  second,  local  cavitation  is  a  small  area  of  low  pressure 
generally  found  to  occur  at  the  fluid- structure  interface.  Both  types  of  cavitation  can 
figure  in  very  prominently  to  the  overall  response  of  the  ship  during  an  UNDEX  event 
and  are  important  factors  that  must  be  accounted  for  in  the  simulation  process  [Ref.  11]. 


1.  Bulk  Cavitation 

In  an  UNDEX  event  the  blast,  and  likewise  the  shock  wave,  propagate  in  a 
spherical  expanding  circle  from  the  detonation  point.  Figure  3  is  a  two-dimensional 
representation  of  the  geometries  involved  in  a  typical  UNDEX  event.  The  incident 
pressure  wave,  a  compressive  pressure  wave,  is  first  to  strike  the  target.  As  this  shock 
wave  reaches  the  free  surface  it  is  reflected  at  the  boundary  as  a  rarefaction  wave,  which 
means  that  the  water  flow  is  directed  opposite  to  the  direction  of  propagation.  This 
tensile  pressure  wave  contributes  to  the  creation  of  bulk  cavitation  due  to  the  marked 
reduction  in  image  pressure  once  the  incident  pressure  wave  has  decayed.  The  pressure 
decay  rate  is  defined  in  accordance  with  Equation  (1).  This  point  is  termed  “cut-off’. 
The  cavitation  pressure  can  be  as  low  as  negative  three  to  four  psi  [Ref.  12],  Figure  4 
illustrates  the  shock  wave  profiles  and  the  “cut-off’  time. 

There  may  also  exist  be  a  bottom  reflection  wave,  caused  by  the  reflection  of  the 
shock  wave  off  of  the  ocean  floor,  though  this  type  of  pressure  wave  is  customarily  of 
less  significance  in  the  UNDEX  event  pertaining  to  surface  ships  since  this  type  of 
pressure  wave  is  heavily  dependent  on  the  ocean  floor  characteristics  and  its  proximity  to 
the  target  [Ref.  9], 
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Figure  3.  Underwater  Explosion  Geometry  [after  Ref.  9] 
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Figure  4.  Shock  Wave  Pressure  Profile  with  Cut-off  Time  [after  Ref.  9] 


The  bulk  cavitation  area  is  formed  due  to  the  water’s  inability  to  support  the 
negative  pressure  resulting  from  the  tensile  forces  of  the  reflected  incident  pressure  wave. 
The  water  vapor  cavity  that  is  created  consists  of  two  separate  boundary  regions,  an 
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upper  boundary  and  a  lower  boundary.  These  boundaries  are  a  function  of  the  size,  type 
and  depth  of  the  charge  that  is  detonated  in  an  UNDEX  event.  An  example  of  this 
dependence  is  shown  in  Figure  5.  This  series  of  subplots,  generated  for  a  100  lb 
Pentaerythritol  Tetranitrate  (PETN)  charge  at  varying  depths,  illustrates  the  difference  in 
resulting  bulk  cavitation  zones.  The  MATLAB®  code  used  to  generate  this  figure  is 
provided  in  APPENDIX  A. 


Bulk  Cavitation  Region  for  Underwater  Explosion:  100  lb  PETN  Charge  at  25ft 
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.  Bulk  Cavitation  Region  for  Underwaler  Explosion:  100  lb  PETN  Charge  at  50ft 
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Bulk  Cavitation  Region  for  Undeiwater  Explosion:  100  lb  PETN  Charge  at  75ft 
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Figure  5.  Bulk  Cavitation  Zone  for  100  lb  PETN  Charge  Detonated  at  Varying  Depths 


In  order  to  find  the  upper  cavitation  zone  boundary,  which  is  defined  as  the  area 
in  which  the  net  pressure  equals  zero,  Equation  (7)  is  used  in  conjunction  with  Equations 
(8)  and  (9). 
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where  x,y  =  the  horizontal  range  and  the  vertical  depth  of  the  point 
7,  =  standoff  distance  from  the  charge  to  the  point 

r2  =  standoff  distance  from  the  image  charge  to  the  point 

C  =  acoustic  velocity  in  the  water 
D  =  charge  depth 

6  =  decay  constant 

Pa  =  atmospheric  pressure 
y  =  weight  density  of  water 

W  =  charge  weight 

Kl,Al,K2,A2=  shock  wave  parameters  (charge  type  dependent) 


In  order  to  determine  the  lower  cavitation  zone  boundary,  the  decay  rates  of  the 
absolute  pressure  and  the  reflected  wave  must  be  equated.  The  formula  for  this 
calculation  is  shown  in  Equation  (10),  using  the  same  variables  as  in  Equations  (7) 
through  (9). 


r2  -  2D 


G{x'y)=~ter 


D  +  y 


+  ^L(Pi+PA+ry)  =  0 

r2 


Pt ,  the  incident  pressure  at  cut-off,  is  given  by, 
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(11) 


P=P  e 

1 1  1  max1^ 


(r2~rl)  j 
09  J 


Bulk  cavitation  will  continue  to  exist  in  the  area  of  the  water  vapor  cavity  until  its 
absolute  pressure  has  risen  above  zero  psi.  The  bulk  cavitation  area,  when  viewed  during 
an  UNDEX  event  occurring  near  to  the  surface,  can  be  witnessed  as  a  white  flattened 
cardioid-like  shape  just  beneath  the  air-water  interface.  Figure  6  is  a  cross  section  view 
representative  of  the  bulk  cavitation  zone  created  by  a  100  lb  PETN  charge  detonated  75 
feet  below  the  free  surface  of  the  water. 


H  V  / 


Figure  6.  Bulk  Cavitation  Zone  Resulting  from  an  Underwater  Explosion 


2.  Local  Cavitation 

Shock  pressure  pulses  occurring  as  a  result  of  an  UNDEX  event  excite  the  ship 
structure  as  they  impinge  on  the  hull,  causing  dynamic  responses.  As  this  fluid-structure 
interaction  occurs,  the  total  pressure  along  the  hull  becomes  negative.  Unable  to  sustain 
the  tension,  the  water  pressure  reduces  to  vapor  pressure  and  cavitation  occurs.  Taylor 
flat  plate  theory  will  be  used  to  describe  how  the  phenomenon  of  local  cavitation  occurs. 
Figure  7  shows  a  Taylor  flat  plate  subjected  to  a  plane  wave. 
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Shock  Wave  Approaching  Shock  Wave  Reflecting 


Figure  7.  Taylor  Plate  Subjected  to  a  Plane  Wave  [after  Ref.  9] 

The  plate  is  subjected  to  an  incident  shock  wave,  Pt(t) ,  which  is  taken  to  be  a 
planar  wave.  As  is  interacts  with  the  plate,  the  reflected  pressure  wave,  P2  (?) ,  is  created 
and  reflected  off  of  the  plate.  The  velocity  of  the  plate  is  defined  as  u(t).  Using 
Newton’s  2nd  Law  the  equation  of  motion  for  a  flat  plate  is  written  as 

m^p-  =  Pfi)  +  P2(t)  (12) 

at 

The  velocities  behind  the  incident  shock  wave  and  the  reflected  shock  wave  are 
defined  as  ux  (?)  and  u2  (?) ,  respectively.  From  this  the  plate  interface  between  the  surface 
of  the  plate  and  the  fluid  is  described  as 

u(t)  =  Mj(?)-m2(?)  (13) 

It  can  be  shown  for  a  one-dimensional  wave,  that  the  that  incident  and  reflected 
shock  wave  pressures  reduce  to, 
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Pl(t)  =  pCu^t) 


(14) 


P2(t )  =  pCu2(t )  (15) 

where  p  =  fluid  density  and  C  =  acoustic  velocity  in  the  medium. 

Hence,  Equations  (1),  (14)  and  (15)  can  be  used  to  formulate  the  solution  of  the 
reflected  pressure  wave  equation  as  follows: 

m=^P-pcu(t)  (i6) 

Equation  (16)  can  then  be  re-written  as  a  first  order  linear  differential  equation, 


(  du\  ) 

m  +pCu(t)  =  2Pmme[s)  (17) 

\dt ) 
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Equation  (19)  then  shows  that  as  P  increases  to  a  large  value,  which  represents  a 
lightweight  plate,  the  total  net  pressure  becomes  negative  at  a  very  early  time.  Thus, 
cavitation  occurs  as  the  vapor  pressure  of  water  is  reached.  The  plate  separates  from  the 
fluid  and  attains  a  maximum  velocity  [Ref.  9], 


D.  FLUID-STRUCTURE  INTERACTION 

The  fluid-structure  interaction  between  the  ship’s  hull  and  the  surrounding  water 
is  primarily  in  the  vertical  direction  as  a  result  of  an  underwater  explosion.  It  has  been 
found  that  the  response  of  the  ship  can  be  approximated  with  some  appreciable  degree  of 
accuracy  by  using  the  Doubly  Asymptotic  Approximation  (DAA)  [Ref.  11].  The 
underlying  equation  used  to  define  the  structural  motion  is  the  discretized  differential 
equation, 


[M]  {x(t)}+[C]  {x(t)}+[K]  |x(t)}={f(t)}  (20) 


where  [ M ] ,  [C]  and  [A']  are  the  symmetric  linear  structural  mass,  damping  and 
stiffness  matrices,  and  {/}  is  the  external  force  vector. 

Equation  (20)  represents  the  dynamic  response  of  the  ship  structure.  It  can  be 
thought  of  as  a  balance  of  all  of  the  forces  acting  upon  the  ship’s  structure.  These  forces 
include  the  inertial  forces,  damping  forces,  internal  forces  and  acoustic  fluid  pressure 
forces  [Ref.  13]. 

In  dealing  with  a  submerged  structure  excited  by  an  acoustic  wave,  the  external 
forcing  function  is  equal  to, 


/  -~GAf(pj  +  ps)  +  fD  (21) 

where  p:  is  the  nodal  pressure  vector  for  the  wetted-surface  fluid  mesh  pertaining 
to  the  incident  wave  and  ps  is  the  nodal  vector  corresponding  to  the  scattering  wave. 
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The  dry-structure  applied  force  vector  is  f D ,  the  transformational  matrix  relating  the 
nodal  surface  forces  is  denoted  as  G ,  while  Af  is  the  diagonal  area  matrix  associated  with 
the  elements  in  the  fluid  mesh  [Ref.  13]. 

The  DAA  is  the  preferred  method  for  solution  of  this  problem  since  it  accounts 
for  the  approximation  of  both  early  time  (high  frequency)  and  late  time  (low  frequency) 
motions  [Ref.  14],  The  First  Order  Doubly  Asymptotic  Approximation  (DAAi)  is  used 
for  long  cylindrical  shell  structures  such  as  a  surface  ship  or  submarine  since  the  DAA  is 
not  valid  in  the  cavitation  region.  The  DDAi  equation  is  expressed  as, 


[Mf  ]  {ps}+pc  [  Af  ]  {ps}=pc  [Mf  ]  {us} 


(22) 


where {us}  is  the  scattered  wave  fluid  particle  velocities  vector,  {ps}  is  the 
scattered  wave  pressure  vector,  [Mf]  is  the  fluid  mass  matrix,  [Af~\  is  the  diagonal  area 
matrix  of  the  fluid  mesh,  p  is  the  fluid  density,  and  c  is  the  sound  velocity  of  the  fluid 
[Ref.  15], 

In  the  early  time  response  Equation  (22)  reduces  to  ps  =  pcus ,  a  plane  wave 
approximation,  while  for  the  late  time  response,  it  reduces  to  Afps  =  Mfus ,  a  virtual  mass 
approximation  [Ref.  16]. 

Since  this  method  allows  for  the  solution  of  the  fluid-structure  interaction  in  terms 
of  wetted-surface  response  only,  {/} ,  the  excitation  of  the  wetted-surface  structure  by 

the  incident  shock  wave  is  given  by  Equation  (23),  where  p,  and  ps  correspond  to  the 
incident  pressure  and  scattered  shock  wave  pressure  vectors,  respectively  [Ref.  17], 

f  =  GAf(pI+ps )  (23) 

The  scattered  wave  fluid  particle  velocities  are  then  tied  to  the  structure  response 
through  the  following  relationship: 
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Gt x  =  Uj  +  us 


(24) 


where  T  represents  matrix  transpose. 

By  means  of  mathematical  manipulation  of  the  aforementioned  equations  the 
resulting  DDA  Interaction  Equations  are  found  to  be: 

[M]  {x}  +  [C]  {x}  +  [K]  {x}  =  -[G]  [A,  ](Pl  +  ps )  (25) 

Wf-l{p,}  +  PWf]{p,}  =  pc{Mt}([GT  {x}-{u,})  (26) 

Equations  (25)  and  (26)  leave  two  unknowns,  x  and  ps ,  which  can  be  solved 
using  a  staggered  solution  scheme  [Ref.  14], 
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II.  SHIP  SHOCK  MODELING  AND  SIMULATION 


A.  SHIP  MODEL 

The  finite  element  model  of  the  USS  WINSTON  S.  CHURCHILL  (DDG-81)  that 
was  used  as  the  basis  of  the  research  reported  in  this  paper  was  built  by  the  DDG-51 
Class  ship  design  firm,  Gibbs  and  Cox,  Inc.  [Ref.  18].  Figure  8  shows  the  Flight  IIA 
model  alongside  the  original  DDG-51  model  used  during  the  DDG-53  simulation  effort. 


Figure  8.  DDG-81  (a)  and  DDG-53  (b)  Finite  Element  Models  [from  Ref.  18] 

The  major  modifications  that  were  made  to  the  Flight  IIA  Arleigh  Burke  Class 
destroyer  initial  production  design  were  as  follow: 

—  Extension  of  the  transom  by  five  feet 

—  Replacement  of  both  Vertical  Launch  System  (VLS)  handling  cranes  by 
six  additional  VLS  cells,  for  a  total  of  96  VLS  cells 

—  Replacement  of  the  5”/54  caliber  gun  with  the  5”/62  caliber  gun 

—  Increased  thickness  of  scantlings  amidships 

—  Installation  of  dual  helicopter  hangers 

Detailed  descriptions  of  the  alterations  made  to  DDG-81  are  provided  in  Figure  9. 
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Figure  9.  Illustration  of  Alterations  Made  to  the  Flight  IIA  Arleigh  Burke  Class 

Destroyer  [from  Ref.  19] 
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For  each  one  of  the  three  shots,  Gibbs  and  Cox,  Inc.  provided  a  separate  finite 
element  model  to  the  Naval  Postgraduate  School  shock  simulation  team.  Details  such  as 
the  liquid  tank  levels,  exact  weapons  load-out,  temporarily  installed  equipment  and  even 
the  number  of  personnel  onboard  at  the  time  of  each  shot,  were  accounted  for  in  order  to 
obtain  the  most  accurate  model  possible  for  simulation  of  the  ship  shock  trials.  The 
complex  finite  element  model  included  many  additional  improvements  over  previous 
models,  such  as  more  realistic  mass  distribution  through  the  use  of  a  significantly 
superior  weight  tape. 

Figure  10  shows  some  of  the  simplified  spring-mass  models  that  were  developed 
from  existing  detailed  equipment  models.  This  was  done  for  items  that  could 
significantly  influence  the  ship  response  at  the  keel,  bulkheads  and  sensor  locations.  The 
gas  turbines,  main  reduction  gears  and  5”/62  gun  are  some  examples  of  critical  items  that 
were  incorporated  into  the  overall  ship  model  [Ref.  18]. 
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Figure  10.  Equipment  Models  in  the  Finite  Element  Model  of  DDG-81  [from  Ref.  18] 
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The  nominal  mesh  size  of  the  finite  element  model  was  27  in  x  48  in.  The  level 
of  detail  and  complexity  of  the  ship  model  are  shown  in  Figure  11,  a  cut-away  view 
looking  from  the  centerline  toward  the  port  side,  astern  of  Frame  300. 


Figure  11.  Cut-away  View  of  the  DDG-81  Finite  Element  Model  [from  Ref.  18] 


Table  1  provides  a  list  of  some  of  the  key  properties  associated  with  the  DDG-81 
finite  element  models. 


Table  1.  List  of  Finite  Element  Model  Properties  [from  Ref.  18] 


Number  of  Nodes 

40,514 

Number  of  Degrees  of  Freedom 

243,084 

Number  of  Beam  Elements 

49,397 

Number  of  Thin  Shell  Elements 

48,662 

Number  of  Lumped  Masses 

92,541 
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Figure  12  gives  a  general  representation  of  where  the  charges  were  located  during 
the  at  sea  ship  shock  trial  performed  on  DDG-81.  The  exact  locations  and  charge  sizes 
will  not  be  discussed  herein  as  these  parameters  are  classified.  Shot  3  was  the  closest  in 
proximity  to  the  ship  when  its  charge  was  detonated.  Inasmuch  as  this  was  the  most 
severe  shot  of  the  three,  the  Shot  3  data  will  be  used  in  most  of  the  comparisons  since  it 
most  closely  approximates  an  actual  surface  ship  UNDEX  event  experienced  during 
combat. 


Figure  12.  DDG-81  Ship  Shock  Trial  Shot  Geometry 
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B.  SHOCK  TRIAL  SIMULATION 

An  overview  of  the  process  used  to  conduct  the  ship  shock  trial  simulations  is 
shown  in  Figure  13.  First  the  finite  element  model  provided  by  Gibbs  and  Cox,  Inc.  was 
converted  from  a  MSC/NASTRAN  input  deck  into  a  nonlinear  dynamic  analysis  code 
(LS-DYNA)  keyword  file.  A  fluid  mesh  was  created  in  TrueGrid,  a  high  quality 
hexahedral  mesh  program  by  XYZ  Scientific  Applications,  Inc.  During  March  through 
June  2003,  Schneider  [Ref.  8]  developed  the  fluid  mesh  that  is  currently  being  employed 
for  all  shock  trial  simulations  at  NPS.  Next  the  fluid  mesh  and  converted  model  were  run 
joined  together  in  LS-DYNA,  which  was  coupled  with  the  Underwater  Shock  Analysis 
code  (USA).  This  is  where  the  actual  shock  simulation  is  conducted.  The  node  output 
data  generated  by  LS-DYNA  was  transferred  to  Ceetron’s  GLview  for  post-processing. 
Finally,  UERD  Tools  was  used  to  compare  the  measured  ship  shock  trial  data  with  the 
shock  simulation  response  results.  A  more  detailed  explanation  of  the  NPS  modeling  and 
simulation  process  follows. 


1.  Pre-Processing 

After  receiving  the  finite  element  model  it  was  translated  into  LS-DYNA 
keyword  format.  A  corresponding  fluid  mesh  for  each  shot’s  model  was  built  in 
TrueGrid.  The  industry  standard  LS-DYNA  software,  which  is  commonly  used  to 
analyze  the  dynamic  response  of  large  structures,  including  those  coupled  to  fluids  was 
chosen  as  a  primary  means  in  which  to  perform  the  simulations.  It  is  a  non-linear  three- 
dimensional  analysis  code  that  performs  the  time  integration  for  the  structure. 


2.  Underwater  Shock  Analysis  Code 

The  USA  code  [Ref.  14]  was  used  to  calculate  the  transient  response  of  the  ship’s 
wetted-surface  structure  to  an  incident  shock  wave.  USA  is  a  boundary  element  code  that 
solves  the  ship’s  structure  interaction  equations  using  the  DAA  formulation  given  in 
Equation  14.  As  previously  stated,  by  using  the  DAA  approach,  the  response  is  modeled 
solely  in  terms  of  the  wetted-surface  variables.  This  eliminates  the  need  for  a  separate 
fluid  volume.  This  technique  has  been  shown  to  work  well  for  a  submerged  structure, 
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such  as  a  submarine,  but  has  some  difficulty  addressing  ship  shock  phenomena  close  to 
the  air- water  interface  due  to  the  addition  of  bulk  cavitation  associated  with  the  surface 
ship  UNDEX  event.  To  overcome  this  problem,  a  finite  element  model  of  the 
surrounding  fluid  elements  was  created  to  properly  account  for  the  presence  of  bulk 
cavitation  within  the  UNDEX  environment  so  that  the  calculations  could  be  performed. 
In  recent  work  completed  by  Hart  [Ref.  20],  it  was  concluded  that  the  surrounding  fluid 
mesh  must  extend  radially  outward  from  the  hull  to  a  radius  equal  to  the  maximum  depth 
of  the  lower  cavitation  boundary.  Accordingly,  the  DAA  boundary  is  truncated  at  the 
outer  surface  of  the  surrounding  fluid  mesh  [Ref.  9]. 


Figure  13.  NPS  Modeling  and  Simulation  Process  Flow  Chart  [from  Ref.  8] 
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The  USA  code  is  formed  from  three  parts,  a)  the  Fluid  Mass  Processor 
(FLUMAS),  b)  the  Augmented  Matrix  Processor  (AUGMAT),  and  c)  the  Time 
Integration  Processor  (TIMINT)  [Ref.  11]. 


a.  FLUMAS 

The  FLUMAS  processor  creates  the  fluid  mass  matrix  for  the  ship’s 
wetted-surface  structure  in  an  infinite,  in  viscid  and  incompressible  fluid.  Based  on  user- 
defines  inputs  which  include:  fluid  mesh  and  element  definitions,  frees  surface  location, 
working  medium  fluid  properties  and  atmospheric  properties,  the  FLUMAS  processor 
calculates  the  number  of  independent  coordinates  required  to  define  the  hull’s  structural 
and  fluid  degrees  of  freedom  (DOF)  on  the  wetted-surface  [Ref.  17,  22],  In  addition,  the 
FLUMAS  processor  generates  the  directional  cosines  for  the  normal  pressure  force  and 
the  nodal  weights  for  the  fluid  element  pressure  forces  [Ref.  16,  21]. 


b.  AUGMAT 

The  AUGMAT  processor  is  where  the  fluid  and  structural  matrices  are 
linked  together.  The  output  from  the  FLUMAS  processor,  specifically  the  symmetric 
fluid  mass  matrix,  is  sent  to  the  AUGMAT  processor  for  use  along  with  the  LS/DYNA 
generated  structural  mass  matrix,  to  create  input  matrices  for  the  TIMINT  processor.  The 
combination  of  these  matrices  within  the  same  file  makes  for  a  more  efficient  manner  in 
which  the  TIMIT  processor  is  able  to  access  the  data  [Ref.  22], 


c.  TIMINT 

The  final  processor  in  the  USA  code  is  the  TIMINT  processor.  It  then 
compiles  the  output  information  from  the  AUGMAT  processor  and  uses  this  data  to 
execute  the  direct  integration  of  Equations  (17)  and  (18).  These  are  the  structural  and 
fluid  interactions,  respectively.  The  TIMINT  processor  solves  the  fluid  equations 
whereas  the  LS-DYNA  processor  solves  the  structural  equations.  Both  of  these  equation 
sets  are  solved  at  every  time  step  by  using  an  unconditionally  stable  staggered  integration 
scheme.  The  TIMINT  processor  output  data  is  saved  as  a  binary  history  file,  (D3THDT), 
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and  as  an  ASCII  file,  (NODOUT).  Thus,  a  time  history  of  displacement,  velocity  and 
wetted-surface  pressure  is  recorded  for  those  nodes  that  were  previously  designated  in  the 
LS-DYNA  keyword  input  file.  Since  the  TIMINT  processor  is  the  most  time  intensive  in 
the  entire  simulation  process,  response  data  information  is  only  retained  for  those  nodes 
that  have  been  chosen  based  on  their  correlation  to  actual  sensor  locations  from  the  ship 
shock  trials  [Ref.  23]. 


3.  Post-Processing 

The  results  obtained  from  the  LS-DYNA  and  USA  codes  are  then  transferred  into 
a  graphical  post-processing  software  package  for  further  conversion  of  the  data  into  a 
visual  representation  of  the  ship  shock  trial  simulation  response  data.  This 
transformation  accomplished  in  Ceetron’s  GLview  Pro  Suite.  The  GLview  output  is  then 
exported  to  the  UERD  Tools  software  where  velocity  time  history  response  plots  are 
generated  for  comparison  of  the  simulation  against  the  measured  ship  shock  trial  data. 


a.  GLview 

Ceetron’s  GLview  Pro  Suite  is  a  commercial  application  that  caters  to  the 
thee-dimensional  visualization  and  interactive  animation  of  simulations  run  on  large 
complex  Finite  Element  models.  GLview  has  the  ability  to  directly  import  binary  and 
ASCII  type  data  files  generated  by  the  LS-DYNA/USA  processors.  Possessing  not  only 
the  capability  of  three-dimensional  model  visualization  but  also  an  ability  to  create  time- 
dependent  data  plots,  GLview  Pro’s  animation  software  is  able  to  display  time-dependent 
results  in  both  scalar  and  vector  formats  for  the  stresses,  strains,  displacement,  velocities 
and  accelerations  within  the  fluid-structure  model  [Ref.  24],  Unfortunately,  GLview  Pro 
is  unable  to  directly  import  ship  shock  trial  data  for  comparison.  Thus  the  ASCII  history 
files  for  each  sensor/node  location  must  be  extracted  from  the  LS-DYNA  NODOUT  file 
and  individually  exported  to  the  UERD  Tools  data  analysis  and  plotting  program. 


29 


b.  UERD 

Underwater  Explosions  Research  Department  (UERD)  is  a  RTD&E 
organization  within  the  Naval  Surface  Warfare  Center,  Carderock  Division.  The  data 
analysis  program,  UERD  Tools,  is  a  custom  software  package  designed  specifically  for 
the  analysis  of  ship  shock  trial  data.  This  extremely  versatile  software  contains  many 
features  for  the  manipulation  and  filtering  of  raw  data  as  well  as  the  conditioning  of 
imported  data  of  various  formats.  UERD  Tools  enables  the  user  to  create  high  quality 
data  plots  of  shock  response.  As  the  final  step  in  the  NPS  modeling  and  simulation 
process,  the  LS-DYNA/USA  simulation  data  is  imported  in  the  ASCII  type  file  format 
generated  in  GLview  Pro.  Prior  to  comparison  of  the  simulation  data  against  the 
measured  sensor  data  the  time  steps  of  all  response  frequency  curves  are  normalized  and 
scaled  to  ensure  proper  fit  of  dimensional  units.  Lastly,  the  actual  data 
analysis/correlation  is  conducted  using  built-in  analysis  tools  within  the  UERD  Tools 
program. 
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IV.  DATA  COLLECTION  AND  ANALYSIS 


The  methods  described  within  this  section  were  used  in  the  data  processing  and 
error  correlation  of  all  of  the  sensor  and  simulation  considered  in  the  series  of  studies 
presented  in  this  paper. 


A.  SHOCK  RESPONSE  DATA  PROCESSING 

In  order  to  properly  compare  the  simulation  data  and  the  measured  data  each  set 
must  be  analyzed  for  anomalies  and  inherent  errors  embedded  within  the  data.  High 
frequency  “noise”  and  low  frequency  “drift”  are  two  such  factors.  They  must  first  be 
minimized  by  established  methods  so  that  their  influence  does  not  skew  the  data 
comparison. 

1.  High  Frequency  “Noise” 

The  sensors  used  in  the  measurement  and  recording  of  actual  ship  shock  trial  data 
not  only  collect  the  desired  frequency  response  but  also  gather  unwanted  high  frequency 
“noise”.  These  additional  frequencies,  which  are  well  beyond  the  interest  range  for 
UNDEX  events,  tend  to  clutter  the  data.  The  unfiltered  data,  shown  in  red  in  Figure  14 
has  a  less  defined  frequency  curve  as  compared  to  the  low-pass  filtered  data,  in  blue,  for 
the  same  sensor.  The  time  history  plot  for  this  velocity  meter,  V2010V  was  taken  from 
the  Shot  2  data  set.  By  using  the  low-pass  filtering  technique,  all  of  the  frequencies 
greater  than  250  Hz  were  removed,  leaving  a  much  cleaner  plot. 

The  aforementioned  process  of  noise  reduction  has  been  an  accepted  practice  for 
some  time,  however  it  has  only  recently  been  postulated  that  the  same  procedure  should 
be  applied  to  the  simulated  data  as  well.  A  statistical  study  based  on  233  accelerometer 
measurements  indicated  that  the  simulation  data,  when  low-pass  filtered  at  250Hz, 
correlated  much  better  with  the  low-pass  filtered  raw  data  for  the  same  sensor  [Ref.  8], 
Table  2  shows  a  summary  of  the  statistical  results  of  this  study.  Consequently,  all 
response  data  comparisons  in  this  study  were  low  pass  filtered  at  250  Hz. 


31 


DDG-81  SHOT  2 


V2010V  -  ICeel  Sensor 


Time  (msec) 

LJn  11  ltered  - Low-Pass  Filtered 


Figure  14.  Comparison  of  Unfiltered  and  Low-Pass  Filtered  Sensor  Data 


Table  2.  Summary  of  Filtered  and  Unfiltered  Simulation  Data  [from  Ref.  8] 


Shock  Trial  Data 
(Filtered) 

Simulation  Data 
(Unfdtered) 

Simulation  Data 
(Filtered) 

Mean 

26.225 

82.985 

34.297 

Variance 

520.229 

5775.711 

606.426 

Standard  Deviation 

22.809 

75.998 

24.626 

2.  Velocity  Response  “Drift” 

For  the  analysis  conducted  in  this  work,  two  types  of  sensor  data  were  used, 
namely  that  collected  from  accelerometers  and  velocity  meters.  Though  they  were 
designed  to  capture  the  transient  response  motion  of  the  system,  these  sensors  are 
routinely  used  to  gather  data  for  up  to  2500  msec  during  an  UNDEX  event. 
Consequently,  the  sensors  acquire  a  larger  range  of  frequencies  than  are  desired. 
Velocity  meters  require  seismic  correction,  an  integration  process,  to  correct  their  error. 
The  drift  associated  with  data  taken  from  accelerometers  is  a  result  of  the  integration 
process  that  transforms  it  into  velocity  response  data. 
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The  time  response  history  for  sensor  A4100V  data  shown  in  Figure  15  was  taken 
from  the  Shot  2  data  set.  A  gradual  trailing  off  of  the  sensor’s  time  history  data  is  seen 
after  the  first  250  msec.  This  trend  away  from  the  zero  equilibrium  point  is  the  drift  in  the 
sensor  data. 


Shot  Number  10622 


Created  03 -J  iin-200  1 


Time  (msec) 

■  A41QOV^ 

Figure  15.  Accelerometer  Output  Data  Prior  to  Drift  Compensation  being  Applied 

There  have  been  numerous  studies  conducted  concerning  the  problem  of  sensor 
drift.  In  this  study,  the  built-in  “Drift  Compensation”  function  within  UERD  Tools  was 
chosen  as  a  means  to  eliminate  the  drift  present  in  the  measured  accelerometer  data. 
Though  this  technique  does  not  always  produce  a  time  response  history  entirely  free  from 
drift,  it  maintains  the  magnitude  of  the  response  and  does  not  introduce  a  phase  shift.  It 
follows  a  set  algorithm  as  opposed  to  some  other  curve  fitting  processes  that  require  the 
skill  and  judgment  of  the  user  to  identify  the  point  where  the  drift  is  introduced  into  the 
record. 

Figure  16  shows  an  example  of  velocity  response  data  acquired  from  same 
accelerometer  A4100V  after  it  has  been  integrated  and  modified  using  the  UERD  Tools 
Drift  Compensation  function. 
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Figure  16.  Accelerometer  Output  Data  After  Drift  Compensation  has  been  Applied 

Table  3  provides  results  from  a  study  comparing  the  vertical  velocity  response 
data  taken  from  accelerometers  versus  that  obtained  from  velocity  meters.  Overall,  there 
was  a  much  better  correlation  between  the  simulated  data  and  the  measured  data  when 
using  the  accelerometer  data.  In  order  to  minimize  error  introduced  into  the  data  by 
physical  drift  of  the  velocity  meters,  accelerometer  data  was  used  whenever  possible. 


Table  3.  Average  Comprehensive  Russell’s  Error  Factor  [from  Ref.  8] 


SENSOR  TYPE 

Average  Comprehensive  Russel’s  Error  Factor 

SHOT  1 

SHOT  2 

SHOT  3 

OVERALL 

ACCELEROMETER 

0.1845 

0.1434 

0.1910 

0.1730 

VELOCITY  METER 

0.2269 

0.2196 

0.2315 

0.2260 
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B. 


DATA  ANALYSIS  AND  COMPARISON 


1.  Sensor  and  Node  Location 

There  were  approximately  620  sensors  installed  in  USS  WINSTON  S. 
CHURCHILL  during  the  summer  of  2001  ship  shock  trials  conducted  at  sea 
approximately  100  nautical  miles  off  the  coast  of  Mayport,  Florida.  Strain  gauges, 
velocity  meters  and  accelerometers  were  installed  in  DDG-81  to  capture  the  time  history 
response  data  during  each  of  the  three  UNDEX  events.  Figures  17  and  18  show  the 
locations  where  the  vertical,  athwartship  and  longitudinal  sensors.  In  order  to  facilitate 
sensor  correlation  during  the  analysis,  specific  nodes  corresponding  to  sensor 
installations  were  built  into  the  structural  model  in  many  instances.  If  a  sensor  location 
did  not  exactly  correspond  to  a  node,  the  closet  node  was  found  from  the  finite  element 
model.  The  selected  nodes  were  then  designated  in  the  LS-DYNA  input  deck  as  nodes 
for  which  to  retain  time  history  response  data  for  comparison. 

Typically  the  vertical  velocity  response  is  analyzed  for  an  UNDEX  event,  since  it 
is  the  principal  response  direction.  Accordingly,  for  the  analysis  of  the  Combat 
Information  Center  (CIC)  area,  the  vertical  velocities  were  chosen  for  comparison. 
However,  in  the  second  part  of  the  analysis  that  is  presented,  the  athwartship  velocities 
were  compared  to  investigate  whether  or  not  the  athwartship  simulations  were  accurate. 


Figure  17.  Sensor  Locations  Depicted  in  Profile  View  of  DDG-81  [from  Ref.  19] 
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Figure  18.  Sensor  Locations  Depicted  in  Top  View  of  DDG-81  [from  Ref.  19] 

2.  Error  Measurements 

Though  always  somewhat  subjective,  there  must  exist  away  to  quantify  how  well 
two  data  sets  correlate  to  one  another.  As  has  been  the  case  in  previous  studies 
concerning  the  simulation  of  the  DDG-81  shock  trials,  this  paper  shall  also  use  Russell’s 
error  factor  [Ref.  25,  26]  as  an  unbiased  measure  of  error  between  the  simulated  data  and 
the  measured  data. 

In  order  to  calculate  the  Russell’s  error,  first,  two  variables  are  defined  as, 


and 


5=Z/2(o2 


(27) 


(28) 


where  fi(i)  and  f2(i)  are  the  measured  and  predicted  response  magnitudes  at  each 
time  step,  which  is  denoted  as  i.  Using  the  variables  A  and  B  from  Equations  (27)  and 
(28),  the  relative  magnitude  error  of  the  correlation  is, 
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(29) 


III  = 


(A-B) 

4ab 


From  Equation  (29)  the  magnitude  error  is  calculated  as, 


RM  =  s/gn(m)log10(l+  |  m  |) 


The  phase  error  is  found  as  follows, 


p=k*k 


(30) 


(31) 


where  0  is  the  normalized  unit  vector  of  the  transient  response.  The  phase 
correlation  between  the  two  data  sets  can  be  computed  as, 


p  = 


C 

4ab 


(32) 


where  C  is  defined  by, 


(0/2(0 


i= 1 


(33) 


The  phase  error  is  calculated  as, 


RP  = 


cos  \p ) 


n 


(34) 
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Equations  (30)  and  (34)  are  used  in  conjunction  with  Equation  (35)  to  determine 
the  comprehensive  error. 


RC  = 


( RM2+RP 2) 


(35) 


Now  that  the  correlation  has  been  defined  in  terms  of  a  comprehensive  error 
factor,  a  range  must  be  set,  delineating  what  will  be  deemed  an  acceptable  span  of  error 
values.  Though  there  is  no  definitive  number  that  characterizes  a  “satisfactory” 
correlation  between  the  data  sets,  the  values  listed  in  Table  4  have  been  used  as  the 
acceptance  criteria  in  both  the  earlier  DDG-53  and  current  DDG-81  ship  shock  trial 
simulation  projects  [Ref.  21]. 


Table  4.  Russell’s  Comprehensive  Error  Factor  Acceptance  Criteria 


RC  <0.15 

EXCELLENT 

0.15  <RC<  0.28 

ACCEPTABLE 

RC  >  0.28 

POOR 

Figure  19  is  a  plot  of  the  data  set  that  was  used  in  determining  the  criteria 
presented  in  Table  4.  Notice  that  in  some  cases  a  comparison  with  a  RC  =  0.25  or  0.26 
was  considered  poor  while  conversely  some  plots  having  correlations  as  high  as  0.33  or 
0.34  were  given  an  acceptable  rating. 
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Figure  19.  Russell’s  Error  Criteria  Determination  Data  [from  Ref.  25] 


The  criteria  established  in  Table  4  were  suggested  to  be  a  valid  measure  of 
acceptance  based  on  comparison  of  500  msec  in  length  using  fully  conditioned  velocity 
response  data  comparisons.  The  data  used  in  these  comparisons  was  subjected  to  drift 
correction  and  low-pass  filtering  at  250  Hz  as  previously  described.  The  acceptance 
criteria  from  Table  4  were  determined  to  be  valid  for  only  the  aforementioned  data 
processing  method. 
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V.  SIMULATION  RESULTS 


During  the  DDG-81  ship  shock  trials  a  total  of  three  underwater  explosions  or 
“shots”  took  place.  All  three  of  these  shots  have  been  simulated  using  the  NPS  modeling 
and  simulation  process.  The  primary  focus  of  this  thesis  is  to  further  validate  the  DDG- 
81  simulation  process  developed  at  NPS  by  investigating  the  shipwide  athwartship 
velocity  response  and  the  localized  vertical  velocity  response  in  the  Combat  Information 
Center  area.  Additionally,  the  effects  of  using  a  new  set  of  proportional  damping 
coefficients  will  be  studied  and  compared  with  those  used  in  the  DDG-53  simulation 
effort  conducted  at  the  NPS. 

A.  CIC  AREA  VERTICAL  VELOCITY  RESPONSE  DATA 

The  Combat  Information  Center  (CIC)  is  one  of  the  focal  points  in  ship’s 
operational  life.  The  main  objective  of  the  CIC  personnel  is  to  ensure  combat  readiness 
by  acting  as  a  central  hub  for  the  gathering,  processing,  and  dissemination  of  all 
Command  and  Control,  Communications,  Computers  and  Intelligence  (C4I)  data  and 
information  throughout  the  ship.  With  all  of  its  electronic  equipment,  communications 
devices  and  weapons  systems  consoles,  determining  the  motion  response  values  for  CIC 
was  a  NAVSEA  priority  during  the  DDG-81  ship  shock  trials  program  [Ref.  19].  The 
addition  of  upgraded  consoles  such  as  the  AN/UYQ-70(V)  Advanced  Display  System 
and  the  physical  rearrangement  of  the  CIC  layout  helped  spur  the  effort  to  quantify  the 
response  that  the  equipment  and  watchstanders  would  be  subjected  to  during  an  UNDEX 
event.  Figure  20  displays  the  layout  of  the  CIC  as  modeled  and  the  location  of  select 
sensors  within  the  compartment.  Figure  21  shows  the  location  of  CIC  with  respect  to  the 
ship’s  profile. 
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Figure  20.  DDG-81  CIC  Console  and  Sensor  Locations  [from  Ref.  18] 


500  450  400  350  300  250  200  150  100  50  0  -50 

FRAME 


Figure  21.  Location  of  Combat  Information  Center  [after  Ref.  27] 
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Vertical  velocity  comparisons  of  the  simulated  response  data  and  the  ship  shock 
trial  data  were  made  for  all  three  of  the  shots.  Table  5  lists  the  sensor  types  and  locations 
of  those  that  were  evaluated  in  this  series  of  comparisons.  All  sensors  analyzed  in  this 
section  were  deck  mounted. 


Table  5.  CIC  Vertical  Velocity  Response  Sensor  Locations 


Sensor 

Compartment 

Grid  ID 

X  (inches) 

Y(inches) 

Z(inches) 

General  Location  Description 

A4005V 

CIC 

212054 

3888 

0 

390 

BASE  OF  CENTER  COMPARTMENT  STANTION 

A4025V 

CIC 

212075 

3744 

54 

390 

FOUNDATION  OF  AAWC  CONSOLE 

A4100V 

CIC 

211924 

3888 

-270 

390 

BTWN  GFCS  &  OPTICAL  SIGHT  CONSOLES 

A4101V 

CIC 

212068 

3936 

27 

MK  124  CONSOLE  BASE 

A4102V 

CIC 

212042 

3936 

-27 

UNDER  MK  124  CONSOLE 

A4104V 

CIC 

212156 

3984 

216 

FOUNDATION  OF  OJ-446C/SLQ-32(V) 

A4106V 

CIC 

212025 

3744 

-54 

FOUNDATION  OF  Q-70  CONSOLE 

A4108V 

CIC 

212153 

3840 

216 

390 

BTWN  RSC  &  TIC  CONSOLES 

A4109V 

CIC 

211924 

3888 

-270 

390 

BTWN  SWS  &  ASUWC  CONSOLES 

A4110V 

CIC 

211973 

3744 

-162 

390 

FOUNDATION  OF  ASWCSO  CONSOLE 

A4111V 

CIC 

212155 

3936 

216 

390 

FOUNDATION  OF  Q70  CONSOLE 

A4408V 

CIC 

211979 

4032 

-162 

390 

BTWN  FOUNDATION  OF  LC01  &  LC02  CONSOLE 

A4409V 

CIC 

211926 

3984 

-270 

390 

FOUNDATION  OF  ATDC-1 

A5503V 

CIC 

211974 

3792 

-162 

390 

AT  L6S,  CENTER  OF  FOUNDATION 

A2104V 

CIC  ANNEX 

222240 

0 

390 

AT  BULKHEAD  174  CENTELINE 

A2101V 

CIC  PROJECTION  RM 

212058 

4080 

0 

390 

AT  BULKHEAD  126  CENTERLINE 

A4105V 

CIC  PROJECTION  RM 

212031 

4032 

-54 

390 

FOUNDATION  OF  CLSD 

A2102A 

CIC  PROJECTION  RM 

212058 

0 

390 

AT  BULKHEAD  126  CENTERLINE 

A2105AI 

CIC  ANNEX 

222240 

3504 

0 

390 

AT  BULKHEAD  174  CENTELINE 

A2106F 

CIC  ANNEX 

222240 

3504 

0 

390 

AT  BULKHEAD  174  CENTELINE 

A2103F 

CIC  PROJECTION  RM 

212058 

4080 

0 

390 

AT  BULKHEAD  126  CENTERLINE 

While  the  true  magnitudes  of  the  simulation  data  comparison  contained  both 
positive  and  negative  values,  indicating  simulated  responses  that  were  both  smaller  and 
larger  than  the  measured  magnitudes  of  the  sensor  response  data,  for  ease  of  plotting,  all 
magnitudes  errors  are  plotted  as  their  absolute  value.  The  true  calculated  magnitudes  are 
found  in  the  corresponding  data  tables  for  each  set  of  plots.  Figure  22  is  a  plot  of  the 
complete  data  set  for  all  three  shots  used  in  the  CIC  area  analysis. 

In  all  but  a  few  exceptions,  the  vertical  velocity  response  values  fall  into  the 
excellent  or  acceptable  range.  Even  those  falling  outside  the  acceptable  region  are  just 
barely  greater  than  the  0.28  cut-off  value,  and  do  not  necessarily  constitute  an  undesirable 
correlation.  The  magnitude  error  is  consistently  low  throughout  the  data  set,  while  it  is 
the  relationship  of  the  phase  that  inevitably  drives  the  overall  comprehensive  error  higher 
in  most  cases. 
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Russell's  Comprehensive  Error  Factor 
CIC  Analysis  (Shots  1,  2  &  3) 


- 0.15 - 0.28  •  Shot  1  •  Shot  2  •  Shot3 


Figure  22.  Russell’s  Comprehensive  Error  Factor  for  CIC 


1.  Shot  1 

Shot  1  was  located  furthest  from  the  ship,  making  it  the  least  severe.  This  shot 
was  detonated  forward  of  the  port  bow,  well  off  of  its  intended  mark,  which  was  abeam 
the  port  side  of  the  ship.  Consequently  it  had  the  most  extreme  aspect  as  was  shown  in 
Figure  12.  Based  on  previous  studies  of  the  DDG-81  ship  shock  trial  simulations  it  was 
surmised  that  this  asymmetric  geometry  would  not  negatively  impact  the  course  of  the 
current  study. 


a.  Error  Comparison 

Using  a  fluid  mesh  model  that  extended  down  to  the  cavitation  depth  of  75 
ft,  an  average  Russell’s  Comprehensive  error  factor  was  previously  found  to  be  0.2162 
during  the  ship  wide  vertical  sensor  analysis  of  Shot  1  [Ref.  8],  In  comparison  the  mean 
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value  of  the  localized  CIC  area  analysis  for  shot  1  was  only  0.1978.  This  improvement 
in  Russell’s  Comprehensive  error  factor  was  anticipated  since  a  general  trend  of 
improvement  as  the  z-direction  coordinate  of  the  sensor  location  increases  with  respect  to 
the  ship’s  keel,  or  baseline  was  previously  discovered  [Ref.  8],  Figure  23  shows  a 
graphical  representation  of  the  Russell’s  Comprehensive  error  factor  while  Table  6 
provides  the  supporting  data. 


Russell's  Comprehensive  Error  Factor 
CIC  Analysis  (Shot  1) 


- 0.15 - 0.28  •  Shot  1 

Figure  23.  Russell’s  Error  Factor  for  CIC  (Shot  1) 
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Table  6.  Russell’s  Error  Factors  for  CIC  (Shot  1) 


Simulation  runtime  =  500  msec  | 

SHOT  1 

Sensor 

NODE 

Mounting 

Type 

Location  (in)* 

738in  Cavitation  -  Medium 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

X 

y 

z 

RM 

RP 

RC 

A4005V 

212054 

Deck 

3888 

0 

390 

0.0876 

0.1637 

0.1645 

A4025V 

212075 

Deck 

3744 

54 

390 

-0.1318 

0.1522 

0.1784 

A4100V 

211924 

3888 

-270 

390 

-0.1168 

0.1566 

0.1732 

A4101V 

212068 

HsSB 

27 

390 

0.0087 

0.1687 

0.1497 

A4102V 

212042 

Deck 

3936 

-27 

390 

0.1321 

0.2234 

0.2300  1 

A4104V 

212156 

Hs^BI 

216 

390 

-0.0018 

0.2525 

0.2238 

A4106V 

212025 

Hs^SI 

3744 

-54 

390 

-0.1906 

0.1903 

0.2387 

A4108V 

212153 

Deck 

3840 

216 

390 

0.1201 

0.3270 

0.3087 

A4109V 

211924 

Hs^Si 

-270 

390 

-0.1297 

0.1685 

0.1885 

A4110V 

211973 

hessb 

3744 

-162 

390 

0.0487 

0.1472 

0.1374 

A4111V 

212155 

Deck 

3936 

216 

390 

-0.0762 

0.2701 

0.2487 

A4408V 

211979 

4032 

-162 

390 

0.2086 

0.2501 

0.2886 

A4409V 

211926 

HEE31! 

3984 

-270 

390 

0.0719 

0.1589 

0.1545 

A2104V 

222240 

Deck 

3504 

0 

390 

-0.0670 

0.1270 

0.1272 

Deck 

4080 

0 

390 

0.0210 

0.1736 

0.1551 

-0.01520 

2.92980 

2.96700 

>0.28 

Poor 

Sum(E(XA2)) 

0.18558 

0.61683 

0.63021 

<0.15 

Excellent 

Mean 

■mmmi  M 

0.19532 

0.19780 

Standard  Deviation  | 

0.11513 

0.05643 

0.05564 

*  Referenced  to  the  G&C  NASTRAN  Model  with  coordinate  origin  located  at  the  stern. 


b.  Velocity  Plot 

The  following  plots  are  of  the  vertical  velocity  comparison  conducted 
between  the  measured  ship  shock  trial  data  and  the  simulation  data  at  typical  sensor 
locations  and  their  corresponding  nodes.  Figure  24  shows  the  time  history  response  of 
the  aft  CIC  bulkhead  which  had  one  of  the  best  overall  correlation  factors  with  an  RC  = 
0.1272  while  Figure  25  shows  the  time  history  response  of  a  sensor,  A4108V,  with  one  of 
the  poorest  correlations  with  an  RC  =  0.3087.  Even  so,  this  correlation  value  is  just 
beyond  the  acceptable  range,  but  still  well  within  one  standard  deviation  of  the  mean. 
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Vertical  Velocity  (ft/sec) 


D  DO-81  SHOT  1  -  CIO 

Grid  22240-vz  (A2 1 04V) 


Bulkhead  174  (x=35()4  y=C)  z=390) 


Time  (msec) 

|  Ship  Shock  Trial  IS]  PS  Simulation  | 

Figure  24.  Deck  Sensor  A2104V 

DDG-81  SHOT  1  -  CIO 

Grid  212153  -vz  (/\4 1 08  V) 


R SC /TIC  Consoles  (x=384()  y=2  1  6  z=39Q) 


Time  (msec) 

Ship  Shock  Trial  - N]  PS  Simulation 


Figure  25.  Deck  Sensor  A4108V 
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2.  Shot  2 

As  shown  in  Figure  12,  Shot  2  was  detonated  directly  abeam  the  starboard  side  of 
DDG-81  during  the  ship  shock  trials.  The  charge  severity  was  the  intermediary  value  of 
the  three  shots,  with  its  standoff  distance  lying  nearer  the  ship  than  it  had  in  Shot  1 . 


a.  Error  Comparison 

In  the  shipwide  comparison  of  vertical  velocity  response  for  Shot  2  the 
mean  Russell’s  Comprehensive  error  factor  was  determined  to  be  0.1912  [Ref.  8].  The 
mean  value  of  the  localized  CIC  area  analysis  for  Shot  2  was  found  to  be  0.2006.  In  this 
case  the  results  are  comparable  with  those  obtained  from  the  shipwide  analysis.  Figure 
26  is  a  graphical  representation  of  the  Russell’s  Comprehensive  error  factor  for  Shot  2. 


Russell's  Comprehensive  Error  Factor 
CIC  Analysis  (Shot  2) 


- 0.15 - 0.28  •  Shot  2 


Figure  26.  Russell’s  Error  Factor  for  CIC  (Shot  2) 


Contrary  to  the  slight  increase  in  overall  RC  value,  the  data  distribution 
throughout  the  shot  has  a  better  accuracy  and  precision  associated  with  it.  Three  of  the 
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sensors  evaluated  fall  in  the  excellent  range  and  only  one  in  the  poor  correlation  category. 
Figure  27  of  sensor  A2104V,  which  had  the  best  correlation  at  RC  =  0.0972,  shows  the 
similitude  of  both  curves  in  magnitude  as  well  as  in  phase.  The  error  in  magnitude  and 
phase  are,  RM  =  0.0623,  RP  =  0.0902,  respectively.  Further  analysis  conducted  for  a  200 
msec  time  history  response  found  that  the  comprehensive  effort  factor  fell  to  a  value  of 
RC  =  0.0799.  Table  7  provides  a  complete  description  of  the  error  factors  for  Shot  2. 


Table  7.  Russell’s  Error  Factors  for  CIC  (Shot  2) 


Simulation  runtime  =  500  msec  | 

SHOT  2 

Sensor 

NODE 

Mounting 

Type 

Location  (in)* 

738in  Cavitation  -  Medium 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

X 

y 

z 

RM 

RP 

RC 

A4005V 

212054 

Deck 

3888 

0 

390 

0.0727 

0.1815 

0.1733 

A4025V 

212075 

Deck 

3744 

54 

390 

-0.0243 

0.1152 

0.1044 

A4100V 

211924 

Deck 

3888 

-270 

390 

-0.1322 

0.1937 

0.2079 

A4101V 

212068 

Deck 

3936 

27 

390 

0.0995 

0.1985 

0.1968 

A4102V 

212042 

Deck 

3936 

-27 

390 

0.0372 

0.2234 

0.2007 

A4104V 

212156 

Deck 

3984 

216 

390 

-0.0637 

0.2739 

0.2492 

A4106V 

212025 

Deck 

3744 

-54 

390 

-0.0973 

0.2091 

0.2044 

A4108V 

212153 

Deck 

3840 

216 

390 

-0.1260 

0.2862 

0.2771  | 

A4109V 

211924 

Deck 

3888 

-270 

390 

-0.1814 

0.1955 

0.2364 

A4110V 

211973 

Deck 

3744 

-162 

390 

0.0442 

0.2132 

0.1930 

A4111V 

212155 

Deck 

3936 

216 

390 

-0.1970 

0.2737 

0.2988 

A4408V 

211979 

Deck 

4032 

-162 

390 

0.1324 

0.2397 

0.2434 

A4409V 

211926 

Deck 

3984 

-270 

390 

-0.0551 

0.2097 

0.1921 

A2104V 

222240 

Deck 

3504 

0 

390 

-0.0623 

0.0902 

0.0972 

A2101V 

212058 

Deck 

4080 

0 

390 

0.0258 

0.1498 

0.1347 

Russell  Error  Correlation 

Sum(E(X)) 

-0.52750 

3.05330 

3.00940 

>0.28 

Poor 

Sum(E(XA2)) 

0.16282 

0.66429 

0.64998 

<0.15 

Excellent 

Mean 

-0.03517 

0.20355 

0.20063 

Standard  Deviation  | 

0.10151 

0.05528 

0.05746 

*  Referenced  to  the  G&C  NASTRAN  Model  with  coordinate  origin  located  at  the  stern. 


b.  Velocity  Plots 

Figures  27  and  28  are  representative  of  the  result  obtained  from  the  Shot  2 
analysis  of  the  vertical  velocity  response  in  CIC  and  are  provided  as  a  sample  of  the 
complete  set  of  time  history  response  plot  found  in  APPENDIX  B.  The  Russell’s 
Comprehensive  error  factors  for  sensors  A2104V  and  A4025V  are,  RC  =  0.0972  and 
RC  =  0.1044,  respectively. 
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Vertical  Velocity  (ft/sec) 


D  DO-81  SHOT  2  -  CIO 

Grid  222240-vz  (A2 1 04V) 


Bulkhead  174  (x=35C>4  y=C)  z=390) 


Time  (msec) 

Ship  Shock  Trial  - NTS  Simulation 


Figure  27.  Deck  Sensor  A2104V 

DDG-81  SHOT  2  -  CIO 

Grid  212075-vz  (A4025V) 


AAWC  Console  (x=3744  y=54  z=390) 


Time  (msec) 

Ship  Shock  Trial  - NTS  Simulation 


Figure  28.  Deck  Sensor  A4025V 
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3.  Shot  3 

Shot  3  was  the  most  severe  of  the  three  shots  conducted  during  the  DDG-8 1  ship 
shock  trials.  During  this  UNDEX  event  the  charge  was  detonated  at  the  closest  point  to 
the  ship  in  order  to  create  a  two-thirds  design  level  blast.  The  charge  was  located  almost 
exactly  amidships  along  the  port  beam,  as  shown  in  Figure  12. 


a.  Error  Comparison 

In  the  shipwide  comparison  of  vertical  velocity  response  for  Shot  3  the 
mean  Russell’s  Comprehensive  error  factor  was  determined  to  be  0.2114  [Ref.  8],  Once 
again  the  CIC  area  specific  mean  value  was  comparable  but  just  slightly  higher,  with  a 
Russell’s  Comprehensive  error  factor  of  0.2238.  In  the  case  of  Shot  3,  there  were  no 
excellent  correlations  however,  as  shown  in  Figure  29,  the  data  falls  in  a  much  more 
accurate  manner.  All  but  one  sensor  was  deemed  to  have  acceptable  correlation  with 
respect  to  the  acceptance  criteria,  with  sensor  A4108V  possessing  a  marginal  value  of  RC 
=  0.2897  as  the  sole  poor  correlation. 


Russell's  Comprehensive  Error  Factor 
CIC  Analysis  (Shot  3) 


- 0.15 - 0.28  •  Shot  3 


Figure  29.  Russell’s  Error  Factor  for  CIC  (Shot  3) 
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Nevertheless,  all  of  the  results  for  Shot  3  are  tightly  clustered  in  the 
acceptable  range  (with  the  exception  of  sensor  A4108V)  nearer  the  ordinate,  indicating  a 
minimal  deviation  in  magnitude  error.  The  greatest  error  in  magnitude  correlation  was  a 
RM  =  0.1320  for  sensor  A2104V.  Based  on  the  summary  of  Russell’s  error  factors 
presented  in  Table  8,  almost  80%  of  the  sensors  have  a  RM  <0.10. 


Table  8.  Russell’s  Error  Factors  for  CIC  (Shot  3) 


Simulation  runtime  =  500  msec 

SHOT  3  i 

Sensor 

NODE 

Mounting 

Type 

Location  (in)* 

738in  Cavitation  -  Coarse 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

X 

y 

z 

RM 

RP 

RC  1 

A4005V 

212054 

Deck 

3888 

0 

390 

0.1250 

0.2521 

0.2494 

A4025V 

212075 

Deck 

3744 

54 

390 

0.0369 

0.2084 

0.1876 

A4100V 

211924 

Deck 

3888 

-270 

390 

-0.0513 

0.2190 

0.1994 

A4101V 

212068 

Deck 

3936 

27 

390 

0.0773 

0.2520 

0.2336 

A4102V 

212042 

Deck 

3936 

-27 

390 

0.1320 

0.2779 

0.2727  | 

A4104V 

212156 

Deck 

3984 

216 

390 

0.0793 

0.2736 

0.2525 

A4106V 

212025 

Deck 

3744 

-54 

390 

-0.0780 

0.2334 

0.2181 

A4108V 

212153 

Deck 

216 

0.1053 

0.3095 

0.2897 

A4109V 

211924 

Deck 

-270 

-0.0653 

0.1927 

0.1803 

A4110V 

211973 

Deck 

3744 

-162 

390 

0.0409 

0.2075 

0.1874  j 

A4111V 

212155 

Deck 

216 

-0.0412 

0.2951 

0.2641 

A4408V 

211979 

Deck 

MliLlJM 

-162 

0.1274 

0.2727 

0.2667  1 

A4409V 

211926 

Deck 

3984 

-270 

390 

0.0343 

0.2223 

0.1993 

A2104V 

222240 

Deck 

0 

-0.0020 

0.1895 

0.1680 

A2101V 

212058 

Deck 

4080 

0 

0.0295 

0.2219 

0.1896 

Russell  Error  Correlation 

Sum(E(X)) 

0.55010 

3.62760 

3.35840 

>0.28 

Poor 

Sum(E(XA2)) 

0.09239 

0.89701 

0.77374 

<0.15 

Excellent 

Mean 

0.03667 

0.24184 

0.22389 

Standard  Deviation 

0.07182 

0.03752 

0.03947 

*  Referenced  to  the  G&C  NASTRAN  Model  with  coordinate  origin  located  at  the  stern. 


b.  Velocity  Plots 

Figures  30  and  31  are  vertical  velocity  plots  chosen  from  analysis  of  Shot 
3  CIC  sensor  locations.  The  Russell’s  Comprehensive  error  factors  for  sensors  A4110V 
and  A4100V  are,  RC  =  0.1874  and  RC  =  0.1994,  respectively. 
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Vertical  Velocity  (ft/sec) 


D  DO-81  SHOT  3  -  CIO 

Grid  21  1973-vz  (A41  lOV) 


ASWCSO  Console  (x=3744  y=-  1  62  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  30.  Deck  Sensor  A4110V 


DDG-81  Shot  3  -  CIO 

Grid  211 924-vz  (A4 1 OOV) 


GFCS  Console  0=3888  y=-27()  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  31.  Deck  Sensor  A4100V 
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4.  Statistical  Analysis  of  CIC  Velocity  Response 

Table  9  presents  statistical  results  from  each  of  the  three  shots.  This  table  has 
been  included  as  an  overview  of  the  data  that  has  been  presented  with  regard  to  the 
vertical  velocity  response  analysis  conducted  on  the  CIC  area.  As  is  shown  in  Table  9 
there  is  a  high  rate  of  correlation  for  all  three  shots. 


Table  9.  Statistical  Data  for  CIC  Response  Analysis  of  Shots  1,  2  &  3 


Russell's  Comprehensive  Error 

Shot  1 

Shot  2 

Shot  3 

RC  <  0.30 

93% 

100% 

100% 

RC  <  0.28 

87% 

87% 

87% 

RC  <  0.25 

87% 

87% 

67% 

RC  <  0.20 

60% 

40% 

47% 

RC  <  0.18 

53% 

27% 

6% 

RC  <  0.15 

20% 

20% 

6% 

Mean  RC 

0.1978 

0.2006 

0.2239 

Standard  Deviation 

0.0564 

0.0576 

0.0395 

Mean  +  Standard  Deviation 

0.2542 

0.2582 

0.2634 

Data  within  One  Standard  Deviation 

87% 

87% 

73% 

The  preceding  results  obtained  from  the  CIC  response  data  indicates  that  the  NPS 
modeling  and  simulation  methodology  does  in  fact  consistently  produce  satisfactory 
results  as  compared  to  the  measured  data. 


B.  SHIPWIDE  ATHWARTSHIP  VELOCITY  RESPONSE  DATA 

The  primary  response  of  an  UNDEX  event  is  in  the  vertical  direction. 
Accordingly,  the  athwartship  direction  and  the  longitudinal  direction  responses  are 
significantly  smaller  in  magnitude  than  those  in  the  vertical  direction.  For  this  reason  the 
vertical  response  of  the  system,  which  is  the  ship  in  this  case,  has  always  been  the  focus 
of  previous  analysis  conducted  in  this  area.  In  the  following  study  the  athwartship 
motion  response  of  the  DDG-81  was  simulated  and  compared  to  the  measured  ship  shock 
trial  data  in  a  similar  manner  to  that  previously  discussed.  The  goal  of  this  study  was  to 
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ascertain  whether  or  not  the  NPS  modeling  and  simulation  methodology  accurately 
captured  the  more  subtle  athwartship  velocity  response  as  well  as  it  did  the  primary 
vertical  velocity  response  of  the  ship.  Data  from  all  three  shots  was  incorporated  into  the 
athwartship  analysis. 

1.  Error  Comparison 

As  a  result  of  its  secondary  nature  in  the  over  response  of  the  ship,  the  athwartship 
motion  of  the  ship  is  not  as  well  documented  as  the  vertical  response.  Of  the  over  600 
sensors  installed  during  the  ship  shock  trial  in  DDG-81  only  about  10%  were  used  to 
collect  athwartship  response  data.  Table  10  is  the  list  of  sensors  used  in  the  athwartship 
analysis.  Though  this  is  a  relatively  small  set  of  data  points,  the  sensor  locations  chosen 
were  well  distributed  throughout  the  ship. 


Table  10.  Athwartship  Velocity  Response  Sensor  Locations 


Sensor 

Compartment 

Grid  ID 

X  (inches) 

Y(inches) 

Z(inches) 

General  Location  Description 

A2001A 

PASSAGE 

120217 

5328 

0 

82 

VERTICAL  CENTERLINE  STIFFNER 

A2102A 

CIC  PROJECTION  RM 

212058 

4080 

0 

390 

AT  BULKHEAD  126  CENTERLINE 

A2110A 

RADAR  ROOM  #1 

414953 

4059 

0 

722.8 

ON  BULKHEAD 

A2105AI 

CIC  ANNEX 

222240 

3504 

0 

390 

AT  BULKHEAD  174  CENTELINE  | 

A2117AI 

RADAR  ROOM  #2 

414367 

3504 

0 

702 

AT  BULKHEAD  174 

A2238AI 

PORT  MAST 

416419 

3504 

135 

848 

MAST  LEG  PORT 

A2241A 

STARBOARD  MAST 

416269 

3504 

-135 

848 

MAST  LEG  STARBOARD 

A2015AI 

AUX  MACH  RM  #2 

230461 

2952 

0 

85.5 
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Sensors  A2001  and  A2021  were  originally  included  in  the  study  but  do  not  appear 
in  the  final  analysis  for  Shot  1  and  Shot  2.  These  two  sensors  correspond  to 
accelerometers  located  in  the  bow  and  at  the  stem  of  the  ship.  They  exhibit  very  poor 
correlation  characteristics  as  compared  to  the  rest  of  the  data.  The  data  points  for  these 
sensors  fell  well  of  the  chart  shown  in  Figure  32,  which  is  a  comparison  of  Russell’s 
Comprehensive  error  factor  for  all  three  shots.  Considered  to  be  outliers,  they  were 
ultimately  excluded. 
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Russell's  Comprehensive  Error  Factor 
Athwartship  Analysis  (Shots  1,  2  &  3) 
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Figure  32.  Russell’s  Comprehensive  Error  Factor  for  Athwartship  Response 


Table  1 1  contains  the  complete  list  of  data  used  for  the  analysis  and  computation 
of  Russell’s  error  factors  for  the  athwartship  velocity  response. 

2.  Velocity  Plots 

The  figures  that  follow  provide  some  examples  of  the  athwartship  velocity  plots. 
A  complete  set  of  athwartship  velocity  plots  are  provided  in  APPENDIX  C.  Figure  33 
shows  a  time  history  response  plot  for  sensor  A2015A  from  Shot  2.  The  correlation  for 
this  accelerometer  located  nearly  amidships  near  the  keel  has  a  value  of  RC  =  0.2847. 
This  correlation  is  just  outside  the  acceptable  range. 


56 


Table  11.  Russell’s  Error  Factor  for  Athwartship  Response 
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DDG-8 1  SHOT  2 

Grid  23046 1  -vy  (A2015AI) 


Keel  (x=2952  y=Q  z=85.5) 


Ship  Shook:  Xrial  Measured  Data  - N]  PS  Simulation 


Figure  33.  Keel  Sensor  A2015AI 

The  next  set  of  figures  illustrates  some  of  the  deficiencies  discovered  while 
analyzing  the  athwartship  response  data.  Analysis  of  the  early  time  response  for  keel 
sensor  A2033A,  shown  in  Figure  34,  revealed  that  it  suffered  from  serious  drift  error 
within  the  first  100  msec.  The  standard  drift  compensation  technique  had  been  applied, 
with  the  result  being  displayed  in  Figure  35.  There  is  a  definite  trailing  off  of  the 
measured  shock  trial  data  curve  away  from  the  abscissa.  Without  manual  manipulation  of 
the  shock  trial  curve  it  was  impossible  to  generate  a  valid  comparison  of  the  simulated 
data.  Similar  in  nature  to  the  drift  correction  issues  experienced  in  the  CIC  area  analysis, 
this  problem  appears  to  be  magnified  by  the  much  smaller  magnitudes  that  are  witnessed 
in  the  athwartship  response.  In  most  cases,  the  magnitudes  in  the  athwartship  direction 
are  on  the  order  of  one  magnitude  smaller  when  compared  to  the  corresponding  vertical 
direction  response  for  the  same  sensor  location. 
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DDG-81  SHOT  1 

Grid  320764-vy  (A2033A) 


Keel  (x  l  152  y=Q  z=196) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data  - NTS  Simulation 


Figure  34.  Keel  Sensor  A2033A 

Shot  Number  10621 


A2033A 


Time  (msec) 

Taw  Sensor  Data  - Drift  Compensated 


Figure  35.  Sensor  A2033A:  Application  of  Drift  Compensation 
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DDG-8 1  SHOT  2 

Grid  41641 9-vy  (A22 3  8AI) 


Bulkhead  (x=3504  y=135  z=848) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  36.  Bulkhead  Sensor  A2238AI 


The  bulkhead  sensor  shown  in  Figure  36,  A2238AI,  appears  to  have  good  initial 
correlation  in  phase,  but  quickly  goes  out  of  phase.  For  the  500  msec  time  history  plot 
shown,  the  RP  =  0.5305,  with  the  RC  =  0.4940.  The  simulated  data  generally  under¬ 
predicted  the  response  after  the  initial  peak  with  the  RM  =  0.1712.  In  general  this  was  a 
poor  correlation,  even  though  both  of  the  curves  appear  to  be  of  almost  same  wave 
pattern. 


3.  Comparison  Results 

Overall  the  results  in  the  athwartship  direction  were  found  to  be  much  less 
promising  than  those  in  vertical  velocity  comparison.  Using  the  same  250Hz  low-pass 
fdtering  and  standard  drift  correction  via  the  UERD  Tools  built  in  function,  the  mean 
correlation  in  the  athwartship  direction  was  determined  to  be  RC  =  0.3922;  well  beyond 
the  RC  =  0.28  acceptable  limit.  Table  12  further  summarizes  the  finding  for  athwartship 
velocity  response. 
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Table  12.  Statistical  Data  for  Athwartship  Response  Analysis  of  Shots  1,  2  &  3 


Russell's  Comprehensive  Error 

Shot  1 

Shot  2 

Shot  3 

RC  <  0.33 

75% 

13% 

9% 

RC  <  0.30 

50% 

13% 

0% 

RC  <  0.28 

25% 

0% 

0% 

RC  <  0.25 

13% 

0% 

0% 

RC  <  0.20 

0% 

0% 

0% 

Mean  RC 

0.3130 

0.4483 

0.4166 

Standard  Deviation 

0.0766 

0.0865 

0.0732 

Mean  +  Standard  Deviation 

0.3896 

0.5348 

0.4898 

Data  within  One  Standard  Deviation 

88% 

88% 

82% 

The  results  obtained  from  the  small  sampling  of  sensors  that  were  studied 
produced  generally  unacceptable  athwartship  response  correlations  based  on  the  criteria 
established  in  Table  4.  Even  so,  the  Russell’s  Magnitude  error  factor  was  equal  to  value 
of  0.20  or  lower,  considerably  lower  in  most  cases.  This  would  indicate  that  the 
simulation  did  in  fact  accurately  capture  the  range  of  the  motion.  However,  as  in  the  CIC 
vertical  velocity  response  comparison  where  it  had  also  been  witnessed,  though  to  a 
lesser  degree,  the  phase  error  dominated  the  error  correlation  due  to  its  pronounced 
excursion  from  the  measured  data.  In  the  athwartship  direction  the  Russell’s  Phase  error 
factor  was  determined  to  lie  in  the  range  of  0.30  or  higher  in  nearly  all  cases. 

One  of  the  possible  contributors  to  the  less  favorable  correlation  in  the 
athwartship  direction  is  the  inherently  smaller  magnitudes  found  in  the  velocity  response 
as  compared  to  those  in  the  vertical  direction.  With  the  ever  present  problem  of  sensor 
drift,  as  previously  shown  in  Figure  35,  the  induced  error  and  method  of  correction 
impact  the  ultimate  curve  comparison  much  more  significantly  in  the  athwartship 
direction  due  to  the  smaller  range  of  motion  in  the  actual  response  for  a  particular  point. 

Additionally,  the  physical  placement  of  the  sensor  in  some  cases  is  suspect.  In 
review  of  the  sensor  installation  descriptions  it  was  discovered  that  some  sensors  were 
mounted  to  the  web  section  of  stiffeners,  equipment  foundations,  longitudinal  bulkheads 
and  other  locations  off  of  the  true  deck.  It  is  postulated  that  the  placement  of  these 
various  sensors  could  impact  the  phase  response  of  the  actual  sensor  during  the  UNDEX 
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event,  skewing  the  recorded  data  by  measuring  the  motion  of  the  lightweight  component 
that  it  is  affixed  to  rather  than  the  ship  itself.  This  type  of  motion,  at  a  presumably  at 
higher  rate  of  oscillation  would  not  be  present  in  the  computer  simulation  of  that  node 
point  in  the  finite  element  model. 
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V.  SHIP  SYSTEM  DAMPING 


Almost  all  of  the  damping  within  a  structure  is  a  result  of  frictional  energy  that  is 
being  dissipated  at  physical  connection  points  such  as  bolted  or  riveted  joints.  However, 
in  a  ship  the  majority  of  connections  are  welded  rather  than  mechanically  joined,  so  there 
is  much  less  energy  dissipation  through  the  welds.  Ships  do  however  provide  a  viable 
means  for  energy  to  escape  the  system.  This  occurs  through  long  cable  runs,  hangers, 
snubbers  and  out  to  the  fluid  surrounding  the  hull  itself  [Ref.  28]. 


A.  PROPORTIONAL  DAMPING  COEFFICIENTS 

A  study  comparing  the  effects  of  ship  system  damping  effects  was  completed  for 
the  DDG-81  [Ref.  28],  Different  proportional  damping  values  were  applied  to  the  LS- 
DYNA  input  deck  and  simulations  were  conducted  for  Shot  2  of  the  DDG-81  ship  shock 
trials.  A  dense  fluid  mesh  model  was  used  for  the  simulation  in  this  comparison.  The 
time  history  plots  of  two  of  these  sets  of  simulations  employing  different  damping 
coefficients  were  compared  against  the  measured  ship  shock  trial  data  in  the  standard 
manner  which  was  outlined  in  the  previous  chapter  of  this  paper. 

Rayleigh  damping,  a  particular  form  of  proportional  damping,  defines  the 
damping  matrix,  [C],  as 


C  =  a[M]  +  j3[K]  (36) 

in  the  general  expression  for  the  structural  equation  of  motion,  Equation  (37). 


[M]{x}  +  [C]{x}+[F]{x}  =  {F} 


(37) 


The  damping  coefficients  a  and  (3  are  constants.  Equation  (36)  can  be  normalized 
using  mass  normalization  such  that 
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(38) 


mT[cm = [2  c®,  w = 

In  a  complex  system  such  as  a  ship,  the  subscript  r,  which  signifies  the  number  of 
modes,  greatly  exceeds  the  two  modes  necessary  to  determine  the  solution  of  the 
equation.  In  the  case  that  the  system  is  over  determined,  the  coefficients  can  be  found 
using  the  measured  data  and  a  least  squares  curve  fitting  method. 

For  each  mode  of  the  ship  response  the  modal  damping  ratio  is  calculated  using 
Equation  (39). 


a 


+  P0)i 


(39) 


1.  NPS  Damping  Values 

A  new  set  of  damping  coefficient  values  was  determined  by  extensive  analysis  of 
measured  data  taken  from  the  DDG-53  ships  shock  trials.  The  ship  was  divided  into  67 
area  groups  for  the  damping  analysis,  which  included  data  from  773  sensors.  For  the 
frequency  spectrum  of  interest,  0  to  250Hz,  both  the  athwartship  and  vertical  response 
were  measured  and  recorded.  A  least  squares  curve  fit,  as  shown  in  Figure  37,  was 
applied  to  each  area  group.  The  area  groups  were  given  weighted  averages  based  on  the 
number  of  modes  used  in  the  curve  fitting  process  necessary  to  find  a  and  p,  which  are 
shown  in  Tables  13  and  14. 


Table  13.  Weighted  Mean  of  a  [from  Ref.  28] 


Athwartship  Direction 

Vertical  Direction 

18.4 

19.2 

Table  14.  Weighted  Mean  of  P  [from  Ref.  28] 


Athwartship  Direction 

Vertical  Direction 

2.82E-06 

2.09E-06 
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Figure  37.  Modal  Damping  Ratio  at  Area  6,  Athwartship  Direction  [from  Ref.  28] 


Thus,  the  new  NPS  Damping  values  are  defined  as  a  =  19.2  and  (3  =  2.09E-6,  in 
the  vertical  direction,  of  which  shall  be  compared  in  this  study.  Similarly,  in  the 
athwartship  direction  the  NPS  Damping  values  were  found  to  be,  a  =  18.4  and  (3  =  2.82E- 
6.  The  great  disparity  in  the  two  damping  coefficients  indicates  that  the  damping  in  the 
system  is  mass  driven. 


2.  DDG-53  Simulation  Damping  Values 

Table  15  gives  the  damping  coefficient  values  that  were  previously  used  in  the 
modeling  and  simulation  effort  of  DDG-53  and  the  early  DDG-81  investigations. 


Table  15.  Damping  Values  from  the  DDG-53  Simulation  Effort  [from  Ref.  29] 


Damping  Value 

a 

P 

4% 

2.64 

4.99E-05 

8% 

4.93 

9.89E-05 
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The  values  listed  in  Table  15  were  found  by  fixing  the  damping  ratio,  at  two 
particular  frequencies,  namely  5Hz  and  250Hz.  The  complete  curves  were  then 
generated  across  the  frequency  spectrum  using  Equation  (39).  Figures  38  and  39 
illustrate  the  various  damping  curves  considered.  The  points  at  5Hz  and  250Hz  indicate 
where  the  4%  and  8%  damping  curves  were  fixed  to  those  particular  values  of  The 
MATLAB®  code  used  to  generate  the  following  plots  is  provided  as  APPENDIX  D. 


RAYLEIGH  DAMPING 


Figure  39.  Proportional  System  Damping  (Logarithmic  Scale) 
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B.  DAMPING  SYSTEM  COMPARISONS 

In  this  study  the  discussion  shall  be  limited  to  a  comparison  between  the  new 
“NPS  Damping  Values”  as  presented  by  Shin  and  Ham  [Ref.  28],  and  the  fixed  4% 
damping  used  in  the  DDG-53  ship  shock  trial  simulation  analysis  and  the  preliminary 
DDG-81  ship  shock  trial  simulation  analysis  conducted  at  NPS. 

1.  Error  Comparison 

The  following  series  of  velocity  response  plots  compares  the  Rayleigh  damping 
coefficients,  a  and  p,  presented  in  Tables  13  and  14  with  coefficients  that  were  used  in 
previous  studies  conducted  on  the  DDG-53  and  DDG-  81,  which  appear  in  Table  15. 

Russell’s  error  factor  was  once  again  chosen  as  a  means  of  comparing  the 
simulated  velocity  response  data  against  the  measured  actual  ship  shock  trial  data.  For 
the  purpose  of  this  study,  an  established  set  of  acceptance  criteria  was  taken  from  the 
values  presented  in  Table  4. 

2.  Velocity  Plots 

The  following  velocity  plots,  which  were  taken  from  the  analysis  conducted  on 
Shot  2,  and  are  typical  of  the  results  discovered  in  this  portion  of  the  study.  Additional 
time  history  response  plots  are  found  in  APPENDIX  E. 

In  these  figures  the  approximate  location  of  each  sensor  is  indicated  on  the  ship 
accompanying  the  time  history  plots  by  a  red  dot.  The  Russell’s  Comprehensive  (RC) 
error  correlation  factor  was  computed  for  each  sensor.  As  the  velocity  response  plot 
comparisons  in  Figure  40  through  Figure  43  show,  there  is  a  noticeably  closer  correlation 
between  the  NPS  damping  values  and  the  ship  shock  trial  data,  than  with  the  fixed  4% 
damping.  The  mean  RC  for  the  4%  Damping  cases  was  0.25  while  in  comparison  when 
the  new  NPS  damping  values  from  Table  13  and  Table  14  were  used,  the  mean  RC  value 
was  only  0.18.  Recalling  that  by  Russell’s  correlation  criteria,  a  value  below  0.15  is 
considered  an  excellent  correlation,  the  simulations  using  the  new  NPS  damping  values 
consistently  showed  better  correlation  and  an  average  reduction  of  approximately  28%  in 
deviation  from  the  recorded  shock  trial  data  versus  those  using  the  fixed  4%  damping. 
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DDG-81  SHOT  2 

Grid  1 42489-vz  (V2002V) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NI  PS  Simulation  (alpha  =  19.2,  beta  =  2. 09  IT- Ob) 
INI  PS  Simulation  (4%  Damping) 


Figure  40.  Vertical  Velocity  Response:  Deck  Sensor  V2002V 


DDG-81  SHOT  2 

Grid  210894-vz  (V2008VI) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

TNI  PS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

TNI  PS  Simulation  (4%  Damping) 


Figure  41.  Vertical  Velocity  Response:  Deck  Sensor  V2008VI 
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DDG-8 1  SHOT  2 

Grid  330769-vz  (V2035V) 
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Ship  Shock  Trial  Measured  Data 
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NPS  Simulation  (4%  Damping) 


Figure  42.  Vertical  Velocity  Response:  Keel  Sensor  V2035V 


DDG-8 1  SHOT  2 

Grid  222436-vz  (V2125V) 


Time  (msec) 

Ship  Shock  Trial 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09 h- Ob) 
NPS  Simulation  (4%  Damping) 


Figure  43.  Vertical  Velocity  Response:  Bulkhead  Sensor  V2125V 
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3.  Damping  Effects  on  Correlation 

Table  16  provides  the  supporting  data  with  Russell’s  error  factors  computed  for 
all  sensors  used  in  this  portion  of  the  system  damping  comparison  study.  Figure  44  is  a 
graphical  representation  of  the  data  presented  in  Table  16.  The  Russell’s  Comprehensive 
error  factor  for  the  simulations  using  the  NPS  damping  values  all  fall  in  the  highly 
acceptable  range,  with  the  exception  of  two  sensors.  These  two  sensors,  as  in  the  earlier 
studies,  correspond  to  sensors  located  at  the  extremities  of  the  ship,  namely  the  bow  and 
stem.  As  before,  there  is  some  hesitation  in  accepting  these  data  points  that  fall  well 
outside  of  the  pattern  of  the  others  within  their  own  grouping.  These  data  points  should 
be  considered  suspect,  but  are  being  included  in  this  portion  of  the  study  for 
completeness. 


Russell's  Comprehensive  Error  Factor 
Shot  2  -  Complete  Data 


- 0.15 - 0.28  ♦  NPS  Damping  •  4%  Damping 

Figure  44.  Russell’s  Error  Factor  for  DDG-81  Shot  2  (Vertical  Velocities) 


70 


Table  16.  Russell’s  Comprehensive  Error  Factor  Correlation  for  DDG-81  Shot  2 
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Referenced  to  the  G&C  NASTRAN  Model  coordinate  origin  located  at  the  stern. 


Figure  44  shows  that  the  results  from  the  simulations  performed  using  the  fixed 
4%  damping  values  were  only  marginally  acceptable  in  most  cases.  The  fixed  4% 
damping  response  data  is  loosely  grouped  whereas  there  is  considerable  improvement  in 
the  accuracy,  especially  in  terms  of  the  magnitude  correlation  in  the  simulations  using  the 
NPS  damping  values.  This  is  demonstrated  by  the  large  grouping  of  NPS  damping  value 
derived  data  points  nearer  the  ordinate. 

4.  Velocity  Meter  Data 

This  section  examines  the  data  collected  from  select  velocity  meters.  In 
comparison  to  the  overall  data,  using  only  the  velocity  meter  data  shows  an  increase  in 
deviation  between  the  simulated  response  and  measured  ship  shock  trial  data.  In  Figure 
45  there  is  a  noticeable  absence  of  data  points  with  excellent  correlation  that  were 
obtained  using  the  NPS  damping  values. 


Russell's  Comprehensive  Error  Factor 
Shot  2  -  Velocity  Meter  Data 


- 0.15 - 0.28  •  NPS  Damping  •  4%  Damping 


Figure  45.  Russell’s  Comprehensive  Error  Factor  for  Shot  2  (Velocity  Meter  Data) 
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The  weaker  correlation  is  also  indicated  through  the  rise  in  the  Russell’s 
Comprehensive  error  factor  as  shown  in  the  mean  values  listed  in  Table  17.  Part  of  the 
reason  for  the  degradation  in  the  overall  correlation  is  the  inclusion  of  the  two  sensors, 
V2000V  and  V2020V,  which  were  located  at  the  bow  and  stem  of  the  ship,  respectively. 
Yet  a  weaker  correlation  at  the  ship’s  extremities  is  inline  with  the  results  obtained  by 
Schneider’s  [Ref.  8]  shipwide  analysis  of  the  vertical  velocities.  That  work  indicated  that 
there  was  a  direct  correlation  between  the  longitudinal  position  of  a  node  within  the  finite 
element  model  and  the  accuracy  of  the  simulated  data  when  compared  to  the 
corresponding  sensor  data.  The  bow  and  stem  areas  consistently  showed  poor  correlation 
of  the  simulated  data  for  all  three  shots.  An  example  of  this  relationship  is  shown  in 
Figure  46. 


Shot  3  Russell's  Error  vs.  Longitudinal  Position 


Location  (x) 

—  High  Density  Mesh  —  Moderate  Density  Mesh  —  Coarse  Mesh 


Figure  46.  Russell’s  Comprehensive  Error  as  a  Function  of  Position  [from  Ref.  8] 

Therefore,  sensors  V2000V  and  V2020V  were  removed  from  further 
consideration. 
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Table  17.  Damping  Comparison  Results  for  DDG-81  Shot  2  (Velocity  Meter) 
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Referenced  to  the  G&C  NASTRAN  Model  coordinate  origin  located  at  the  stern. 


Considering  only  the  velocity  meter  data,  it  was  found  that  by  using  the  new  NPS 
damping  values,  the  average  improvement  was  approximately  25%  over  the  comparisons 
made  using  the  4%  fixed  damping.  These  results  are  shown  in  Table  18. 


Table  18.  Relative  Improvement  Using  NPS  Damping  Values  (Shot  2) 


Relavtive  Percentage  of  Improvement  in 

Russell's  Comprehensive  Error  Correlation 

Shot  2  (Vertical  Direction) 

Sensor 

Node 

Ship  Compartment 
Location 

Percent 

Relavtive 

Change 

V2002V 

142489 

4th  Deck,  Bow 

32.13% 

V2007V 

210430 

4th  Deck,  Forward 

18.73% 

V2009VI 

210808 

4th  Deck,  Forward 

36.25% 

V2008VI 

210894 

4th  Deck,  Forward 

21.86% 

V2108V 

212196 

1st  Deck,  Bow 

4.32% 

V2010V 

220589 

4th  Deck,  Forward 

15.85% 

V2012VI 

221102 

4th  Deck,  Forward 

12.91% 

V2011VI 

221188 

4th  Deck,  Forward 

25.59% 

V2013V 

221601 

4th  Deck,  Amidships 

-0.26% 

V2124V 

222060 

1st  Deck,  Forward 

-2.25% 

V2125V 

222436 

1st  Deck,  Forward 

24.30% 

V2014V 

230461 

4th  Deck,  Amidships 

14.25% 

V2016V 

242399 

4th  Deck,  Amidships 

20.32% 

V2026V 

312302 

4th  Deck,  Amidships 

14.35% 

V2034V 

330759 

3rd  Deck,  Aft 

-27.43% 

V2035V 

330769 

1st  Deck,  Bow 

31.57% 

V2019V 

340167 

3rd  Deck,  Aft 

17.59% 

Average  Improvement  in  Correlation 

24.58% 

From  investigation  of  the  velocity  meter  data  it  was  clear  that  the  magnitude 
correlation  was  acceptable  in  most  cases,  and  that  the  phase  error  accounted  for  the 
majority  of  the  error.  Thus,  further  analysis  was  conducted  using  a  time  history  of  only 
250  msec,  in  the  order  to  try  and  isolate  the  early  time  response.  Table  19  provides  the 
Russell’s  Comprehensive  error  factors  for  the  shorter  250  msec  comparisons.  Table  20 
shows  the  relative  change  in  Russell’s  Comprehensive  error  factor.  On  average  the  RC 
value  improved  about  30%  when  the  response  time  was  limited  to  the  early  time  frame. 
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Table  19.  Selected  Russell’s  Error  Factors  for  Shot  2  (250  msec) 


Ship  Shock  Simulation  with  Shot  2  Geometry,  Dense  Mesh  and  738  in  Cavitation  Depth 


Sensor 

Node 

Mounting 

Type 

Location  (in)* 

Shock  Trial  Data  vs.  4%  Damping 

Shock  Trial  Data  vs.  NPS  Damping 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

X 

y 

z 

RM 

RP 

RC 

RM 

RP 

RC 

120217 

■USX 

0 

82 

-0.1565 

0.2703 

0.2768 

-0.2817 

0.2599 

0.3397 

■gPESEl 

EBCBI 

4080 

-176 

171 

0.1373 

0.2600 

MiMiLLIiM 

0.1868 

0.1747 

V2013V 

0 

280 

0.1254 

0.1142 

0.1503 

-0.0700 

0.0934 

0.1034 

V2124V 

mzzm 

Bulkhead 

0.0511 

0.2024 

0.1850 

-0.1418 

0.1720 

0.1976 

V2034V 

330759 

Keel 

1152 

-135 

193 

0.0376 

0.1439 

0.1318 

-0.1138 

0.1456 

0.1638 

V2020V 

350052 

Keel 

288 

0 

211 

-0.0962 

0.2610 

0.2456 

-0.2161 

0.2423 

0.2877 

Russell  Error  Correlation 

Sum(E(X)) 

0.09870 

1.25130 

1 .24950 

-0.88640 

1.10000 

1.26690 

>0.28 

Poor 

Sum(E(XA2)) 

0.07235 

0.28324 

0.27872 

0.16798 

0.22066 

0.30526 

<0.15 

Excellent 

Mean 

0.01645 

0.20855 

0.20825 

-0.14773 

0.18333 

0.21115 

Runtime  =  250  msec 

Standard  Deviation 

0.11893 

0.06675 

0.06085 

0.08606 

0.06163 

0.08689 

*  Referenced  to  the  G&C  NASTRAN  Model  coordinate  origin  located  at  the  stern. 


Table  20.  Russell’s  Comprehensive  Error  Factor  for  Early  Time  Response 


Sensor 

Node 

Mounting 

Type 

Location  (in) 

Relative  Change  for 

250  msec  vs.  500  msec 

X 

y 

z 

4%  Damping 

NPS  Damping 

V2000V 

120217 

Keel 

5328 

0 

82 

6.99% 

10.61% 

V2009VI 

210808 

Bulkhead 

4080 

-176 

171 

23.26% 

19.12% 

V2013V 

221601 

Bulkhead 

3504 

0 

280 

3.72% 

33.93% 

V2124V 

222060 

Bulkhead 

3504 

-375 

390 

22.98% 

19.54% 

V2034V 

330759 

Keel 

1152 

-135 

193 

19.83% 

21.81% 

V2020V 

350052 

Keel 

288 

0 

211 

13.06% 

10.62% 

Mean  Relative  Change 

29.63% 

30.97% 

5.  Accelerometer  Data 

When  looking  at  only  the  accelerometer  data  for  the  500  msec  time  history 
response  comparisons,  it  was  determined  that  there  was  an  improvement  31%  over  the 
complete  data  set  and  a  39%  improvement  over  the  velocity  meter  only  data  comparisons. 
The  excellent  precision  in  the  data  points  obtained  from  the  NPS  damping  value 
simulations  is  clearly  shown  in  Figure  47.  Nearly  all  the  points  have  an  excellent 
correlation  rating,  while  those  from  the  4%  fixed  damping  simulations  are  only 
marginally  acceptable  or  have  poor  correlation. 
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Russell's  Comprehensive  Error  Factor  for  Vertical  Velocity  Response 
Shot  2  -  Accelerometer  Data 


- 0.15 - 0.28  •  NPS  Damping  •  4%  Damping 


Figure  47.  Russell’s  Comprehensive  Error  Factor  for  Shot  2  (Accelerometer  Data) 


Figures  48 
from  the  DDG-81 
coefficient  values. 


through  55  show  the  excellent  correlation  of  the  accelerometer  data 
ship  shock  trial  simulations  performed  using  the  new  NPS  damping 
The  Russell’s  error  factors  for  these  figures  are  provided  in  Table  21. 


Table  21.  Damping  Comparison  Results  for  DDG-81  Shot  2  (Accelerometer) 


Ship  Shock  Simulation  with  Shot  2  Geometry,  Dense  Mesh  and  738  in  Cavitation  Depth 


Sensor 

Node 

Mounting 

Type 

Location  (in)* 

Shock  Trial  Data  vs.  4%  Damping 

Shock  Trial  Data  vs.  NPS  Damping 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

SHIP  SHOCK  TRIAL  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

LS-DYNA/USA  DATA  (<250HZ) 

X 

y 

z 

RM 

RP 

RC 

RM 

RP 

RC 

A8516V 

210993 

Bulkhead 

1591 

36 

453 

0.1945 

0.1927 

0.2426 

0.0311 

0.1623 

0.1464 

A2104V 

222240 

Bulkhead 

3504 

0 

390 

0.1535 

0.1079 

0.1663 

-0.0388 

0.0953 

0.0912 

A3565V 

231696 

Bulkhead 

2952 

-81 

317 

0.1817 

0.2328 

0.2617 

0.0290 

0.1675 

0.1506 

A2310V 

320746 

Bulkhead 

1536 

0 

177 

0.1948 

0.2776 

0.3005 

0.0103 

0.2047 

0.1816 

A2116V 

414367 

Deck 

3504 

0 

702 

0.1844 

0.1288 

0.1993 

0.0070 

0.1003 

0.0891 

A2109V 

414953 

Bulkhead 

4059 

0 

723 

0.2963 

0.2395 

0.3376 

0.1170 

0.1647 

0.1790 

A2240V 

416269 

Mast 

3504 

-135 

848 

0.1459 

0.1571 

0.1900 

-0.0626 

0.1335 

0.1307 

A2237V 

416419 

Mast 

3504 

135 

848 

0.1815 

0.1405 

0.2034 

0.0047 

0.0931 

0.0826 

Russell  Error  Correlation 

Sum(E(X)) 

1.53260 

1 .47690 

1.90140 

0.09770 

1.12140 

1.05120 

>0.28 

Poor 

Sum(E(XA2)) 

0.30838 

0.29840 

0.47646 

0.02110 

0.16906 

0.14929 

<0.15 

Excellent 

Mean 

0.19158 

0.18461 

0.23768 

0.01221 

0.14018 

0.13140 

Runtime  =  500  msec 

Standard  Deviation 

0.04594 

0.06065 

0.05922 

0.05333 

0.04117 

0.03994 

*  Referenced  to  the  G&C  NASTRAN  Model  coordinate  origin  located  at  the  stern. 
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Vertical  Velocity  (ft/sec)  Vertical  Velocity  (ft/sec) 


DDG-81  SHOT  2 

Grid  210993-vz  (A8516V) 


I  bulkhead  (x=1591  y=36  z=453) 


Time  (msec) 

Ship  Shook  Xrial  Measured  Data 

TN  PS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
IN  PS  Simulation  (4%  Damping) 


Figure  48.  Bulkhead  Sensor  A8516V 


DDG-81  SHOT  2 

Grid  222240-vz  (A2  1 04V) 


Bulkhead  (x=3504  y=C)  z=390) 


Time  (msec) 

Ship  Shock  Xrial  IV1  earned  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
INI  PS  Simulation  (4%  Damping) 


Figure  49.  Bulkhead  Sensor  A2104V 
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Vertical  Velocity  (ft/sec)  Vertical  Velocity  (ft/sec) 


DDG-81  SHOT  2 

Grid  231696-vz  (A3565V) 


Bulkhead  (x=2952  y=-8  1  z=3  1  7) 


Time  (msec) 

Ship  Shook  Xrial  Measured  Data 

NJ  PS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
N]  PS  Simulation  (4 °/o  Damping) 


Figure  50.  Bulkhead  Sensor  A3565V 

DDG-81  SHOT  2 

Grid  4 1 4367-vz  (A2 1  1 6V) 


Deck  (x=3504  y=0  z=702) 


Time  (msec) 

Ship  Shock  Xrial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

NPS  Simulation  (4%  Damping) 


Figure  51.  Bulkhead  Sensor  Bulkhead 
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Vertical  Velocity  (ft/sec) 


DDG-81  SHOT  2 

Grid  320746-vz  (A2310V) 


Bulkhead  (x=1536  y=C)  z=177) 


Time  (msec) 


Ship  Shook  Xrial  Measured  Data 

NJ  PS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
N]  PS  Simulation  (4 °/o  Damping) 


Figure  52.  Deck  Sensor  A2116V 

DDG-81  SHOT  2 

Grid  414953-vz  (A2109V) 


Bulkhead  (x=4()59  y=Q  z=723) 


Time  (msec) 

Ship  Shock  Xrial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

NPS  Simulation  (4%  Damping) 


Figure  53.  Bulkhead  Sensor  A2109V 
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Vertical  Velocity  (ft/s) 


DDG-81  SHOT  2 

Grid  4 1 6269-vz  (A2240V) 


Mast  (  x=35()4  y=-135  z=848) 


Time  (msec) 

Ship  Shook  Xrial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NPS  Simulation  (4 °/o  Damping) 


Figure  54.  Mast  Sensor  A2240V 

DDG-81  SHOT  2 

Grid  41641 9-vz  (A2237V) 


Mast  (x=35()4  y=135  z=848) 


Time  (msec) 

Ship  Shook  Xrial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

NPS  Simulation  (4%  Damping) 


Figure  55.  Mast  Sensor  A2237V 
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VII.  CONCLUSION  AND  RECOMMENDATIONS 


A.  SUMMARY  OF  FINDINGS 

Based  upon  the  finding  presented  herein,  and  those  referenced  throughout  this 
paper,  it  becomes  evident  that  the  modeling  and  simulation  methodology  developed  at  the 
Naval  Postgraduate  School  Shock  and  Vibration  Computational  Laboratory  is  valid  for 
the  simulation  of  naval  surface  ship  system  response  resulting  from  the  standard  LFT&E 
UNDEX  event  designed  to  test  the  ship  up  to  a  two-thirds  design  limit.  This  thesis 
investigated  the  primary  system  response  in  the  Combat  Information  Center  area  as  well 
as  the  shipwide  secondary  response  in  the  athwartship  direction.  In  addition,  the  vertical 
velocity  response  data  obtained  through  the  use  of  the  new  Rayleigh  damping  coefficients 
developed  at  the  Naval  Postgraduate  School  were  compared  against  existing  damping 
coefficient  values  previously  employed  in  the  DDG-53  ship  shock  trial  simulations.  With 
the  exception  of  the  athwartship  damping  analysis,  which  determined  the  simulation  data 
to  be  only  marginally  acceptable,  the  results  produces  from  this  series  of  parametric 
studies  were  all  highly  favorable  in  nature. 

The  Combat  Information  Center  is  a  critical  compartment  within  a  naval  surface 
combatant.  For  this  reason  it  is  imperative  that  the  watchstanders  and  equipment  located 
in  this  space  be  able  to  not  only  withstand  the  initial  shock  but  also  be  able  to  continue  to 
properly  function  after  being  subjected  to  and  UNDEX  event.  The  nodal  simulation 
conducted  for  the  CIC  sensors  showed  a  solid  correlation  with  the  measured  data  from 
the  DDG-81  ship  shock  trial.  The  results  from  this  localized  area,  deep  in  the  heart  of  the 
ship,  proved  to  be  consistent  with  the  primary  velocity  response  correlations  performed  at 
the  ship  wide  level. 

The  athwartship  velocity  response  results  were  less  promising  than  were 
anticipate,  but  not  truly  unexpected.  The  very  nature  of  the  data  is  much  different  from 
that  of  the  vertical  response.  The  magnitudes  witnessed  in  the  athwartship  direction  are 
routinely  on  the  order  of  a  magnitude  smaller  than  those  in  the  primary  response 
direction.  Although  the  correlation  was  less  acceptable  by  the  Russell’s  Comprehensive 
error  factor  standards  of  RC  <  0.28,  there  appeared  to  be  a  bias,  rather  than  a  random  type 
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error  associated  with  the  excursion  from  the  measured  response  data.  Suspect  in  the 
analysis  was  the  drift  error,  which  seemed  to  have  overwhelmed  the  response  data 
comprised  of  much  smaller  velocities.  Additionally,  the  pronounced  phase  error 
observed  in  the  simulated  athwartship  data  appeared  to  distort  the  overall  correlation. 
Further  investigation  is  highly  recommended  in  this  area,  specifically  with  respect  the 
late  time  phase  error. 

Finally,  this  study  also  concluded  that  the  new  set  of  system  damping  coefficients 
developed  at  NPS  by  Shin  and  Ham  show  a  distinct  improvement  in  correlation  over  the 
previously  employed  4%  fixed  damping.  The  results  of  this  study  support  the  further  use 
of  these  new  Rayleigh  damping  coefficients  in  the  DDG-81  ship  shock  trial  modeling  and 
simulation  effort.  In  doing  so,  it  was  also  confirmed  that  accelerometer  data  was  better 
suited  for  comparison  of  the  simulated  data  against  the  measured  shock  trial  data. 


B.  FUTURE  WORK 

There  is  still  much  work  to  be  done  in  the  modeling  and  simulation  of  UNDEX 
events.  A  few  possibilities  for  further  course  of  study  are  offered  to  that  end.  Having 
validated  the  NPS  methodology  used  in  simulating  the  current  two-thirds  design  limit 
shock  trial,  scaling  the  charge  shock  factor  to  full  scale  shot  at  the  design  limit  or  higher 
is  the  next  logical  step  in  this  simulation  effort.  Successful  completion  of  a  full-scale 
shock  trial  simulation  would  allow  ship  designers  to  glean  valuable  information 
concerning  the  limiting  design  case,  previously  unobtainable  by  conventional  testing 
means.  The  gains  attained  through  the  global  analysis  of  the  DDG-81  ship  shock  trial 
simulation  effort  suggest  focused  study  of  localized  phenomena  experienced  during  an 
UNDEX  event  such  as  whipping.  In  order  to  further  enhance  the  simulation  of  DDG-81 
ship  shock  trials  the  effects  derived  from  the  introduction  of  elasto-plastic  material 
properties  within  the  finite  element  model  of  the  ship  should  be  investigated.  Lastly,  as 
was  previously  mentioned,  the  secondary  system  response  in  the  athwartship  direction 
requires  further  investigation.  Ultimately,  it  is  desired  to  be  able  to  apply  the  modeling 
and  simulation  techniques  that  have  been  developed  here  at  the  Naval  Postgraduate 
School  to  other  ship  classes  through  a  set  of  design  parameters  based  on  the  findings 

ascertained  through  the  investigation  of  the  DDG-53  and  DDG-81  ship  shock  trials. 
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APPENDIX  A.  BULK  CAVITATION  ZONE  PROGRAM 


The  following  program  code  was  written  using  MATLAB®  6.1  release  12.1.  The 

purpose  of  this  program  is  to  compute  the  bulk  cavitation  zone  boundaries  and  create  a 

visualization  of  the  bulk  cavitation  zone. 

%  Computation  of  Bulk  Cavitation  Zone  for  Underwater  Explosions 
%  LT  Jarema  M.  Didoszak,  USN 
%  APR  2003,  Last  Modified  DEC  2003 
%  SVCL,  Naval  Postgraduate  School,  Monterey,  CA 

%  This  program  is  used  to  develop  the  bulk  cavitation 
%  envelope  for  an  underwater  charge  of  PETN  at  various 
%  charge  weights  and  depths. 

clear;  clc; 

%  Parameter  definitions  are  for  PETN  type  charge 

K1  =  24589;  %  Pmax 

A1  =  1.194;  %Pmax 

K2  =  0.052;  %  Decay  constant 

A2  =  -0.257;  %  Decay  constant 

%  Constants 
Pa  =  14.7; 

Gamma  =  62.5/144; 

C  =  5; 

counter  =  0; 

i=i; 

for  W  =  [100,200,300];  %  Equivalent  charge  weights  (lb) 

for  Dl=  [25,50,75];  %  Charge  location  depths  (ft) 

counter  =  counter+1; 

A=zeros(50,1000); 
fory  =  1:51; 
for  x  =  1:1001; 

R1  =  sqrt((Dl  -  (y-l))A2  +  (x-l)A2);  %  Distance  from  charge  to  desired  location  (ft) 

R2  =  sqrt((Dl  +  (y-l))A2  +  (x-l)A2);  %  Distance  from  image  charge  to  desired  location  (ft) 

theta  =  K2  *  ( WA(  1/3))*  ((( WA(  1  / 3  ))/R  1 ) A( A2)) ;  %  Decay  Constant  (msec) 

Pi  =  (Kl*(WA(l/3)/Rl)A(Al))*(exp(-(R2  -Rl)/(C*theta)));  %  Incident  Pressure  Wave  (psi) 

Ph  =  Gamma*(y-1);  %  Hydrostatic  Pressure  at  y  (psi) 

Pr  =  (K1  * (( WA(  1  /3 )/R2) A( A 1 ))) ;  %  Refelcted  Pressure  Wave  (psi) 


F  =  Pi  +  Pa  +  Ph  -  Pr;  %  Upper  Bulk  Caviataion  Boundary 

G 1  =  -Pi/(C  *  theta)  *  ( 1  +(((R2-2  *  D 1  *  ((D 1  +(y- 1  ))/R2))/R  1 )  *  ( A2  *R2/R  1  -  A2- 1 ))) ; 

G2  =  -(Al*Pi/RlA2)*(R2-2*Dl*((Dl+(y-l))/R2)); 


%  Atmpospheric  pressure  (psi) 

%  Weight  density  of  water  (lb/ftA3) 

%  Acoustic  velocity  of  water  (ft/msec) 
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G3  =  Gamma*((D  1  +(y- 1  ))/R2) ; 
G4  =  (Al/R2)*(Pi+Pa  +  Ph); 

G  =  G1  +  G2  +  G3  +  G4; 


%  Lower  Bulk  Cavitation  Boundary 
%  Combine  Bulk  Cavitation  Boundaries 


ifF>  0.001; 
if  G  <  0; 

A(y,x)=  1; 
end 
end 

if  G  >  0; 

A(y,x)  =  1; 
end 
end 
end 

temp(:,:,  counter)  =  A; 
end 

charge=num2  s  tr(  W) ; 

figure(i)  %  Plots  for  100  lb  charge  PETN 
orient  landscape 
hold  on 
subplot(3,l,l) 
spy(temp(:,:,l)) 

title([’Bulk  Cavitation  Region  for  Underwater  Explosion:  ‘,  charge  , 

‘  lb  PETN  Charge  at  25ft’]) 
xlabel( ‘Radius  (ft)’) 
ylabel( ‘Depth  (ft)’) 
subplot(3,l,2) 
spy(temp(:,:,2)) 

title([’Bulk  Cavitation  Region  for  Underwater  Explosion:  ‘,  charge, ... 

‘  lb  PETN  Charge  at  50ft’]) 
xlabel(‘Radius  (ft)’) 
ylabel(‘Depth  (ft)’) 
subplot(3,l,3) 
spy(temp(:,:,3)) 

title([’Bulk  Cavitation  Region  for  Underwater  Explosion:  ‘,  charge, ... 

‘  lb  PETN  Charge  at  75ft’]) 
xlabel( ‘Radius  (ft)’) 
ylabel( ‘Depth  (ft)’) 

i=i+i; 

end 
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APPENDIX  B.  CIC  VERTICAL  VELOCITY  RESPONSE  PLOTS 


A.  SHOT  1 

The  following  velocity  plots  are  from  the  analysis  of  the  Combat  Information 
Center  for  Shot  1.  The  Russell’s  error  factors  for  the  corresponding  data  correlation 
follow  each  figure  caption. 


DDG-81  SHOT  1  -  CIC 

Grid  22240-vz  (A2 1 04  V) 


Bulkhead  174  (x=3504  y=0  z=390) 


Time  (msec) 


Ship  Shook  Trial  - N]  PS  Simulation 


Figure  56.  Deck  Sensor  A2104V:  (RM  =  0.0670,  RP  =  0.1270,  RC  =  0.1272) 
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DDG-8 1  SHOT  1  -  CIC 

Grid  212058-vz  (A2101V) 


Bulkhead  126  ( x^()8()  y=C)  z=390) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  57.  Deck  Sensor  A2101V:  (RM  =  0.0210,  RP  =  0.1736,  RC  =  0.1551) 

DDG-81  SHOT  1  -  CIC 

Grid  2  1 2054-vz  (A4005V) 


Center  of  Compartment  (  x=3888  y=C)  z=390) 


Time  (msec) 

Ship  Shock  Trial  - NTS  Simulation 


Figure  58.  Deck  Sensor  A4005V:  (RM  =  0.0876,  RP  =  0.1637,  RC  =  0.1645) 
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DDG-81  SHOT  1  -  CIC 

Grid  212075-vz  (A4025V) 


AAWC  Console  (x=3744  y=54  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  59.  Deck  Sensor  A4025V:  (RM  =  0.1318,  RP  =  0.1522,  RC  =  0.1784) 

DDG-81  SHOT  1  -  CIC 

Grid  211 924-vz  (A4 1 OOV) 


GFCS  Console  (X=3888  y^270  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  60.  Deck  Sensor  A4100V:  (RM  =  0.1168,  RP  =  0.1566,  RC  =  0.1732) 
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DDG-81  SHOT  1  -  CIC 

Grid  2  1 2068-vz  (A4 1 0 1 V) 


IVIK  1  24  Console  (x=3936  y=27  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  61.  Deck  Sensor  A4101V:  (RM  =  0.0087,  RP  =  0.1687,  RC  =  0.1497) 

DDG-81  SHOT  1  -  CIC 

Grid  2  1 2042-vz  (A4 1 02  \T) 


MK124  Console  LIP  (x=3936  y=27  z=39(l) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  62.  Deck  Sensor  A4102V:  (RM  =  0.1321,  RP  =  0.2234,  RC  =  0.2300) 
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DDG-81  SHOT  1  -  CIC 

Grid  21215  6-vz  (A4 1 04  V) 


SLQ-32(V)  (x=3984  y=2  1  6  z=39C>) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  63.  Deck  Sensor  A4104V:  (RM  =  0.0018,  RP  =  0.2525,  RC  =  0.2238) 

DDG-81  SHOT  1  -  CIC 

Grid  212025-vz  (A4106V) 


Q-70  Console  L2S  (x=3744  y=-54  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  64.  Deck  Sensor  A4106V:  (RM  =  0.1906,  RP  =  0.1903,  RC  =  0.2387) 
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DDG-81  SHOT  1  -  CIC 

Grid  212153  -vz  (A4 1 08  V) 


RSC/TIC  Consoles  (x=384()  y=2  1  6  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  65.  Deck  Sensor  A4108V:  (RM  =  0.1201,  RP  =  0.3270,  RC  =  0.3087) 

DDG-81  SHOT  1  -  CIC 

Grid  211 924-vz  (A4 1 09  V) 


SWS/ASUWC  Consoles  (X=3888  y^270  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  66.  Deck  Sensor  A4109V:  (RM  =  0.1297,  RP  =  0.1685,  RC  =  0.1885) 
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DDG-81  SHOT  1  -  CIC 

Grid  21  1973-vz  (A4  1  1  OV) 


ASWCSO  Console  (x=3744  y=-  1  62  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  67.  Deck  Sensor  A4110V:  (RM  =  0.0487,  RP  =  0.1472,  RC  =  0.1374) 

DDG-81  SHOT  1  -  CIC 

Grid  212  155-vz  (A4 1  1  IV) 


Q-70  Console  LSP  (x=3936  y=216  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  68.  Deck  Sensor  A4111V:  (RM  =  0.0762,  RP  =  0.2701,  RC  =  0.2487) 
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DDG-81  SHOT  1  -  CIC 

Grid  211 979-vz  (A4408V) 


LC01/LC02  Console  (x=4032  y=-  1  62  z  =390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  69.  Deck  Sensor  A4408V:  (RM  =  0.2086,  RP  =  0.2501,  RC  =  0.2886) 

DDG-81  SHOT  1  -  CIC 

Grid  211 926-vz  (A4409V) 


ATDC- 1  (x=3984  y^2VO  x=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  70.  Deck  Sensor  A4409V:  (RM  =  0.0719,  RP  =  0.1589,  RC  =  0.1545) 
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B.  SHOT  2 


The  following  velocity  plots  are  from  the  analysis  of  the  Combat  Information 
Center  for  Shot  2.  The  Russell’s  error  factors  for  the  corresponding  data  correlation 
follow  each  figure  caption. 


DDG-81  SHOT  2  -  CIC 

Grid  222240-vz  (A2 1 04V) 


Bulkhead  174  (x=3504  y=0  z=39(l) 


Figure  71.  Deck  Sensor  A2104V:  (RM  =  0.0623,  RP  =  0.0902,  RC  =  0.0972) 
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Vertical  Velocity  (ft/sec) 


D  DO-81  SHOT  2  -  CIO 

Grid  212058-vz  (A21Q1V) 


Bulkhead  126  (  x=4()8()  y=C)  z=390) 


Time  (msec) 


Ship  Shock  Trial  - NTS  Simulation 


Figure  72.  Deck  Sensor  A2101V:  (RM  =  0.0258,  RP  =  0.1498,  RC  =  0.1347) 

DDG-81  SHOT  2  -  CIO 

Grid  212054-vz  (A4005V) 


Center  of  Compartment  (  x=3888  y=()  z=390) 


Time  (msec) 

Ship  Shock  Trial  - NTS  Simulation 


Figure  73.  Deck  Sensor  A4005V:  (RM  =  0.0727,  RP  =  0.1815,  RC  =  0.1733) 
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DDG-81  SHOT  2  -  CIC 

Grid  212075-vz  (A4025V) 


AAWC  Console  (x=3744  y=54  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  74.  Deck  Sensor  A4025V:  (RM  =  0.0243,  RP  =  0.1152,  RC  =  0.1044) 

DDG-81  SHOT  2  -  CIC 

Grid  211 924-vz  (/\4 1 OOV) 


GFCS  Console  (x=3888  y=-27Q  z=39Q) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  75.  Deck  Sensor  A4100V:  (RM  =  0.1322,  RP  =  0.1937,  RC  =  0.2079) 
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DDG-81  SHOT  2  -  CIC 

Grid  2  1 2068-vz  (A4 1 0 1 V) 


IVIK  1  24  Console  (x=3936  y=27  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  76.  Deck  Sensor  A4101V:  (RM  =  0.0995,  RP  =  0.1985,  RC  =  0.1968) 

DDG-81  SHOT  2  -  CIC 

Grid  212042-vz  (A4102V) 


MKT 24  Console  LIP  (x=3936  y=27  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NPS  Simulation 


Figure  77.  Deck  Sensor  A4102V:  (RM  =  0.0372,  RP  =  0.2234,  RC  =  0.2007) 
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DDG-81  SHOT  2  -  CIC 

Grid  21215  6-vz  (A4 1 04 V) 


SLQ-32(V)  (x=3984  y=216  z=39C>) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  78.  Deck  Sensor  A4104V:  (RM  =  0.0637,  RP  =  0.2739,  RC  =  0.2492) 

DDG-81  SHOT  2  -  CIC 

Grid  212025-vz  (A4106V) 


Q-70  Console  L2S  <x=3744  y=-54  z=39()) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  79.  Deck  Sensor  A4106V:  (RM  =  0.0727,  RP  =  0.1815,  RC  =  0.1733) 
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DDG-81  SHOT  2  -  CIC 

Grid  212025-vz  (A4106V) 


RSC/TIC  Consoles  (x=3840  y=216  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  80.  Deck  Sensor  A4108V:  (RM  =  0.1260,  RP  =  0.2862,  RC  =  0.2771) 

DDG-81  SHOT  2  -  CIC 

Grid  211 924-vz  (A4 1 09  V) 


SWS/ASUWC  Consoles  (x=3888  y^270  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  81.  Deck  Sensor  A4109V:  (RM  =  0.1814,  RP  =  0.1995,  RC  =  0.2364) 
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DDG-81  SHOT  2  -  CIC 

Grid  21  1973-vz  (A4  1  1  OV) 


ASW/CSO  Console  (x=3744  y=-162  z=39C>) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  82.  Deck  Sensor  A4110V:  (RM  =  0.0442,  RP  =  0.2132,  RC  =  0.1930) 

DDG-81  SHOT  2  -  CIC 

Grid  212  155-vz  (A4 1  1  IV) 


Q-70  Console  LSP  (x=3936  y=216  z=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  83.  Deck  Sensor  A4111V:  (RM  =  0.1970,  RP  =  0.2737,  RC  =  0.2988) 
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DDG-81  SHOT  2  -  CIC 

Grid  211 979-vz  (A4408V) 


LC01/LC02  Console  (x^()32  y=-162  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - NTS  Simulation 


Figure  84.  Deck  Sensor  A4408V:  (RM  =  0.1342,  RP  =  0.2397,  RC  =  0.2434) 

DDG-81  SHOT  2  -  CIC 

Grid  211 926-vz  (A4409V) 


ATDC- 1  (x=3984  y^2VO  x=390) 


Time  (msec) 

Ship  Shook:  Trial  - NTS  Simulation 


Figure  85.  Deck  Sensor  A4409V:  (RM  =  0.0551,  RP  =  0.2097,  RC  =  0.1921) 
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C.  SHOT  3 

The  following  velocity  plots  are  from  the  analysis  of  the  Combat  Information 
Center  for  Shot  3.  The  Russell’s  error  factors  for  the  corresponding  data  correlation 
follow  each  figure  caption. 


DDG-8 1  SHOT  3  -  CIC 

Grid  212058-vz  (A2101V) 


Bulkhead  126  (  x=4080  y=0  z=390) 


Time  (msec) 


Ship  Shock  Trial  - NTS  Simulation 


Figure  86.  Bulkhead  Sensor  A2101V:  (RM  =  0.0295,  RP  =  0.2219,  RC  =  0.1896) 
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DDG-8 1  SHOT  3  -  CIC 

Grid  222240-vz  (A2  1 04V) 


Bulkhead  174  (x=35()4  y— O  z=390) 


Time  (msec) 


Ship  Shock  Trial  - NTS  Simulation 


Figure  87.  Deck  Sensor  A2104V:  (RM  =  0.0020,  RP  =  0.1895,  RC  =  0.1680) 


DDG-81  SHOT  3  -  CIC 

Grid  212054-vz  (A4005V) 


Center  of  Compartment  (x=3888  y  O  z=39Q) 


Time  (msec) 


Ship  Shock  Trial  - Nl  PS  Simulated 


Figure  88.  Deck  Sensor  A4005V:  (RM  =  0.1250,  RP  =  0.2521,  RC  =  0.2494) 
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DDG-81  SHOT  3  -  CIC 

Grid  212075-vz  (A4025V) 


AAWC  Console  (x— 3744  y=54  z=39Q) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  89.  Deck  Sensor  A4025V:  (RM  =  0.0369,  RP  =  0.2048,  RC  =  0.1876) 


DDG-81  Shot  3  -  CIC 

Grid  211 924-vz  (A4 1 OOV) 


GFCS  Console  (x=3888  y^270  z=390) 


Time  (msec) 

Ship  Shook  Trial  - NTS  Simulation 


Figure  90.  Deck  Sensor  A4100V:  (RM  =  0.0513,  RP  =  0.2190,  RC  =  0.1994) 
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DDG-81  SHOT  3  -  CIC 

Grid  212068  (A4101V) 


IVIK  1  24  Console  LIP  (x=3936  y =27  z=390) 


Time  (msec) 

Ship  Shook  Trial  - NPS  Simulation 


Figure  91.  Deck  Sensor  A4101V:  (RM  =  0.0773,  RP  =  0.2520,  RC  =  0.2336) 

DDG-81  SHOT  3  -  CIC 

Grid  2  1 2042-vz  (A4 1 02  V) 


MK  124  Console  LIS  (x=3936  y^27  z— 390) 


Time  (msec) 

Ship  Shook  Trial  - NPS  Simulation 


Figure  92.  Deck  Sensor  A4102V:  (RM  =  0.1320,  RP  =  0.2779,  RC  =  0.2727) 
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DDG-81  SHOT  3  -  CIC 

Grid  21215  6-vz  (A4 1 04  V) 


SLQ-32(V)  (x=3984  y=2  1  6  z=390) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  93.  Deck  Sensor  A4104V:  (RM  =  0.0793,  RP  =  0.2736,  RC  =  0.2525) 

DDG-81  SHOT  3  -  CIC 

Grid  212025-vz  (A4106V) 


Q-70  Console  L2S  <x=3744  -y=-54  z=390) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  94.  Deck  Sensor  A4106V:  (RM  =  0.0780,  RP  =  0.2334,  RC  =  0.2181) 
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DDG-81  SHOT  3  -  CIC 

Grid  212153  -vz  (A4 1 08  V) 


RSC/TIC  Consoles  (x=384()  y=216  z=39Q) 


Time  (msec) 

Ship  Shook  Trial  - NTS  Simulation 


Figure  95.  Deck  Sensor  A4108V:  (RM  =  0.1053,  RP  =  0.3095,  RC  =  0.2897) 

DDG-81  SHOT  3  -  CIC 

Grid  211 924-vz  (A4 1 09  V) 


SWS/ASUWC  Consoles  (x=3888  y=-270  z=390) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  96.  Deck  Sensor  A4109V:  (RM  =  0.0653,  RP  =  0.1927,  RC  =  0.1803) 
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DDG-81  SHOT  3  -  CIC 

Grid  21  1973-vz  (A4  1  lOV) 


ASWCSO  Console  (x=3744  y=-162  z=390) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  97.  Deck  Sensor  A4110V:  (RM  =  0.0409,  RP  =  0.2075,  RC  =  0.1874) 

DDG-81  SHOT  3  -  CIC 

Grid  212  155-vz  (A4 1  1  IV) 


Q-70  Console  LHP  (x=3936  y=216  z— 390) 


Time  (msec) 

Ship  Shook  Trial  - NTS  Simulation 


Figure  98.  Deck  Sensor  A4111V:  (RM  =  0.0412,  RP  =  0.2951,  RC  =  0.2641) 
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DDG-81  SHOT  3  -  CIC 

Grid  211 979-vz  (A4408V) 


LC01/LC02  Consoles  (x^()32  y=-  1  <32  z=390) 


Time  (msec) 

Ship  Shook  Trial  - NTS  Simulation 


Figure  99.  Deck  Sensor  A4408V:  (RM  =  0.1274,  RP  =  0.2727,  RC  =  0.2667) 

DDG-81  SHOT  3  -  CIC 

Grid  211 926-vz  (A4409V) 


ATDC- 1  (x=3984  y^270  z=39(l) 


Time  (msec) 


Ship  Shook  Trial  - NTS  Simulation 


Figure  100.  Deck  Sensor  A4409V:  (RM  =  0.0343,  RP  =  0.2223,  RC  =  0.1993) 
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APPENDIX  C.  ATHWARTSHIP  VELOCITY  PLOTS 


A.  SHOT  1 

The  following  velocity  plots  are  taken  from  the  Athwartship  analysis  of  Shot  1 . 
The  Russell’s  error  factors  for  the  corresponding  data  correlation  follow  each  figure 
caption. 


DDG-81  SHOT  1 

Grid  23046 1  -vy  (^A2015^AI) 


Keel  (x=2952  y  O  z=85.5) 


Time  (msec) 


Ship  Shook  Trial  Measured  Data  - NPS  Simulation 


Figure  101.  Keel  Sensor  A2015AI:  (RM  =  0.2403,  RP  =  0.3245,  RC  =  0.3578) 
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DDG-81  SHOT  1 

Grid  320764-vy  (A2033A) 

Keel  (x=l  152  y=Q  z=196) 


Time  (msec) 


■  Ship  Shock  Trial  Measured  Data 


■  NTS  Simulation 


Figure  102.  Keel  Sensor  A2033A:  (RM  =  0.0553,  RP  =  0.5207,  RC  =  0.4641) 


DDG-81  SHOT  1 

Grid  212058-vy  (2A21022A) 


Bulkhead  (x=4080  y=G  z=390) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data  - NTS  Simulation 


Figure  103.  Bulkhead  Sensor  A2102AI:  (RM  =  0.0708,  RP  =  0.2367,  RC  =  0.2189) 
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DDG-81  SHOT  1 

Grid  4 1 6269-vy  (A224 1 2A) 


Bulkhead  (x=3504  y=-135  z=848) 


Ship  Shook  Trial  Measured  Data  - NTS  Simulation 


Figure  104.  Bulkhead  Sensor  A2241A:  (RM  =  0.0847,  RP  =  0.3119,  RC  =  0.2864) 


DDG-81  SHOT  1 

Grid  41641 9-vy  (A223 8AI) 


Bulkhead  (x=3504  y=  1  35  z=848) 


Ship  Shock  Trial  Measured  Data  - NTS  Simulation 


Figure  105.  Bulkhead  Sensor  A2238AI:  (RM  =  0.0179,  RP  =  0.3542,  RC  =  0.3143) 
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DDG-81  SHOT  1 

Grid  414953 -vy  (A2 1  1 OA) 

Bulkhead  (x=-4Q59  y— O  z=723) 


Time  (msec) 


Ship  Shook  Trial  Measured  Data 


■  NTS  Simulation 


Figure  106.  Bulkhead  Sensor  A2110A:  (RM  =  0.1393,  RP  =  0.3381,  RC  =  0.3241) 

DDG-81  SHOT  1 

Grid  222240-vy  (A2 1 05AI) 


Bulkhead  (x=3504  y — O  z=390) 


Ship  Shock  Trial  Measured  Data  - NTS  Simulation 


Figure  107.  Bulkhead  Sensor  A2105AI:  (RM  =  0.0470,  RP  =  0.2873,  RC  =  0.2580) 
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DDG-8 1  SHOT  1 

Grid  4 1 4367-vy  (A2 1  1 7AI) 


Bulkhead  (x=3504  y— O  z=702) 


Time  (msec) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  108.  Bulkhead  Sensor  A2117AI:  (RM  =  0.0311,  RP  =  0.2997,  RC  =  0.2670) 


B.  SHOT  2 

The  following  velocity  plots  are  taken  from  the  Athwartship  analysis  of  Shot  2. 
The  Russell’s  error  factors  for  the  corresponding  data  correlation  follow  each  figure 
caption. 
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DDG-81  SHOT  2 

Grid  414953 -vy  ( A2 1  1 OA) 


Bulkhead  lx-4059  y=Q  z=V23) 


Figure  109.  Bulkhead  Senor  A2110A:  (RM  =  0.1565,  RP  =  0.4387,  RC  =  0.4127) 

DDG-81  SHOT  2 

Grid  212058-vy  (A2102A) 


Bulkhead  (x=4080  y=C)  z=39Q) 


Figure  110.  Bulkhead  Sensor  A2102A:  (RM  =  0.2431,  RP  =  0.3542,  RC  =  0.3807) 
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DDG-8 1  SHOT  2 

Grid  222240-vy  (A2 1 Q5AI) 


Bulkhead  (x=3504  y— O  z=39Q) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  111.  Bulkhead  Sensor  A2105AI:  (RM  =  0.1694,  RP  =  0.4986,  RC  =  0.4667) 

DDG-81  SHOT  2 

Grid  4143 67-vy  (2X2  \  \  72VE) 


Bulkhead  (x=3504  y — O  z=702) 


Ship  Shock  Xrial  Measured  Data  - NPS  Simulation 


Figure  112.  Bulkhead  Sensor  A2117AI:  (RM  =  0.1695,  RP  =  0.5077,  RC  =  0.4744) 


117 


DDG-8 1  SHOT  2 

Grid  41641 9-vy  (A223 8AI) 


Bulkhead  (x=3504  y=135  z=848) 


Time  (msec) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  113.  Bulkhead  Sensor  A2238AI:  (RM  =  0.1712,  RP  =  0.5305,  RC  =  0.4940) 

DDG-81  SHOT  2 

Grid  4 1 6269-vy  (2V22412X) 


Bulkhead  (x=3504  y=-  1  35  z=848) 


Time  (msec) 


Ship  Shock  Xrial  Measured  Data  - NPS  Simulation 


Figure  114.  Bulkhead  Sensor  A2241A:  (RM  =  0.1735,  RP  =  0.5570,  RC  =  0.5170) 
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DDG-8 1  SHOT  2 

Grid  23046 1  -vy  (A2015AI) 


Keel  (x=2952  y=Q  z=8 5.5) 


Ship  Shook  Trial  Measured  Data  - NTS  Simulation 


Figure  115.  Keel  Sensor  A2015AI:  (RM  =  0.1488,  RP  =  0.2847,  RC  =  0.2847) 


C.  SHOT  3 

The  following  velocity  plots  are  taken  from  the  Athwartship  analysis  of  Shot  3. 
The  Russell’s  error  factors  for  the  corresponding  data  correlation  follow  each  figure 
caption. 


119 


DDG-81  SHOT  3 

Grid  12021 7-vy  (A200 1 A) 


Kleel  (x=5328  y— O  z=82) 


Ship  Shook  Trial  Measured  Data  - NTS  Simulation 


Figure  116.  Keel  Sensor  A2001A:  (RM  =  0.4587,  RP  =  0.4203,  RC  =  0.5513) 

DDG-81  SHOT  3 

Grid  23046 1  -vy  (A2015AI) 


Keel  (x=2952  y=0  z=85.5) 


Time  (msec) 


Ship  Shock  Trial  Measured  Dat  - NTS  Simulation 


Figure  117.  Keel  Sensor  A2015AI:  (RM  =  0.2428,  RP  =  0.3212,  RC  =  0.3568) 
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DDG-81  SHOT  3 

Grid  3350052-vy  (A2021A) 


Keel  (x=288  y=0  z=21 1) 


Figure  118.  Keel  Sensor  A2021 A:  (RM  =  0.4226,  RP  =  0.4295,  RC  =  0.5340) 


DDG-81  SHOT  3 

Grid  3307 64-vy  (A2033A) 


Keel  (x=l  152  y— O  z=196) 


Time  (msec) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  119.  Keel  Sensor  A2033A:  (RM  =  0.2102,  RP  =  0.3771,  RC  =  0.3826) 
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DDG-81  SHOT  3 

Grid  212058 -vy  (A2102A) 


Bulkhead  (x=-4080  y— O  z=39Q) 


Ship  Shook  Trial  Measured  Data  - NTS  Simulation 


Figure  120.  Bulkhead  Sensor  A2102A:  (RM  =  0.2265,  RP  =  0.3223,  RC  =  0.3491) 

DDG-81  SHOT  3 

Grid  222240-vy  (A2 1 05AI) 


Bulkhead  (x=3504  y — O  z=390) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data  - NTS  Simulation 


Figure  121.  Bulkhead  Sensor  A2105AI:  (RM  =  0.2691,  RP  =  0.3725,  RC  =  0.4072) 
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DDG-81  SHOT  3 

Grid  414953 -vz  (A2 1  1 OA) 


Bulkhead  (x=^1059  y— O  z— 723) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  122.  Bulkhead  Sensor  A2110A:  (RM  =  0.2156,  RP  =  0.3889,  RC  =  0.3941) 

DDG-81  SHOT  3 

Grid  4143 67-vy  (  A2  \  \  vy^i) 

Bulkhead  (x=3504  y — O  z=720) 


•  Ship  Shock  Xrial  Measured  Data 


■  NPS  Simulation 


Figure  123.  Bulkhead  Sensor  A2117AI:  (RM  =  0.3084,  RP  =  0.3551,  RC  =  0.4168) 
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DDG-8 1  SHOT  3 

Grid  41641 9-vy  (A223 8AI) 


Mast  (x=3504  y=135  z=848) 


Ship  Shook  Trial  Measured  Data  - NTS  Simulation 


Figure  124.  Mast  Sensor  A2117AI:  (RM  =  0.3248,  RP  =  0.3883,  RC  =  0.4487) 

DDG-81  SHOT  3 

Grid  4 1 6269-vy  (A224 1 2A) 


Mast  (x=3504y^!35  z=848) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data  - NTS  Simulation 


Figure  125.  Mast  Sensor  A2241A:  (RM  =  0.2954,  RP  =  0.3741,  RC  =  0.4224) 
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A. thwarts h i p  Velocity  (ft/s) 


DDG-81  SHOT  3 

Grid  320746-vy  (A23  1  1  A) 


Bulkhead  (x-1536  y— O  z=177) 


Ship  Shook  Xrial  Measured  Data  - NPS  Simulation 


Figure  126.  Bulkhead  Sensor  A2311A:  (RM  =  0.1937,  RP  =  0.3040,  RC  =  0.3194) 
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APPENDIX  D.  RAYLEIGH  DAMPING  PROGRAM 


The  following  program  code  was  written  using  MATLAB®  6.1  release  12.1.  The 
purpose  of  this  program  is  to  compute  the  Rayleigh  damping  coefficients  for  the  damping 
comparison  portion  of  this  study. 

%  DDG-81  DAMPING  COEFFICIENTS  COMPARISON 
%  Ralyeigh  Damping  [C]  =  alpha  [M]  +  beta  [K] 

%  FT  Jarema  M.  Didoszak,  USN 

%  JUN  2003,  Modified  OCT  2003  from  3%  to  4%  Damping 
%  SVCF,  Naval  Postgraduate  School,  Monterey,  CA 

clc 

format  long 

%  NPS  Damping 
clear  all 

inputs=XFSREAD('nps,); 
zetal=  inputs  (:, 2); 
omegal=  inputs(:,l)*(2*pi); 

Al=  [  1. /(2* omega  1),  omegal/2]; 

alphabeta  1  =((conj  ( A 1 ')  *  A 1 ) A- 1 )  *  conj  ( A 1 ')  *zeta  1 

W=l:1885; 

regen  1  =alphabeta  1(1 )  ./(2  *  W)+alphabeta  1  (2)  *  W/2 ; 

%  4%  and  4%  damping 
inputs=XFSREAD('44damping'); 
zeta44=  inputs(:,2); 
omega44=  inputs(:,l)*(2*pi); 

A44=  [l./(2*omega44),  omega44/2]; 
alphabeta44=((conj  ( A44’)  *  A44) A- 1 )  *  conj  ( A44’)  *  zeta44 
W=l:1885; 

regen44=alphabeta44(l)./(2*W)+alphabeta44(2)*W/2; 

%  8%  and  8%  damping 
inputs=XFSREAD('88damping’); 
zeta88=  inputs(:,2); 
omega88=  inputs(:,l)*(2*pi); 

A88=  [l./(2*omega88),  omega88/2]; 

alphabeta88=((conj(A88')*A88)A-l)*conj(A88’)*zeta88 

W=l:1885; 

regen88=alphabeta88(l)./(2*W)+alphabeta88(2)*W/2; 
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figure(l) 

plot(W/(2  *pi), regen  1  ,'b') 
hold  on 

plot(W/(2*pi),regen44,'g') 
plot(  W/(2  *pi),regen8  8  ,'r') 

%legend('NPS  DAMPING 
%  ’4%  DAMPING 

%  ’8%  DAMPING  ') 

plot(omega  1 1(2  *pi),zeta  1  ,’b') 
plot(omega44/(2*pi),zeta44,'k.') 
plot(omega8  8/(2  *pi)  ,zeta8  8  ,'k. ') 
axis([0  300  0.0001  0.35]) 
title('RAYLEIGH  DAMPING’) 

xlabel(’FREQUENCY  (Hz)’),  ylabel('DAMPING  RATIO’) 
hold  off 

figure  (2) 

loglog(W/(2*pi),regenl  ,'b') 
hold  on 

loglog(W/(2*pi),regen44,'g') 

loglog(W/(2*pi),regen88,'r') 

%legend('NPS  DAMPING 
%  ’4%  DAMPING 

%  ’8%  DAMPING  ’) 

loglog(omega  1  /(2  *pi)  ,zeta  1  ,'b') 
loglog(omega44/(2*pi),zeta44,'k.') 
loglog(omega8  8/(2  *pi)  ,zeta8  8  ,'k. ') 
axis([2e-l  3e2  5e-3  lei]) 
title(’RAYLEIGH  DAMPING’) 

xlabel(’FREQUENCY  (Hz)’),  ylabel('DAMPING  RATIO’) 
hold  off 
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APPENDIX  E.  DAMPING  VERTICAL  VELOCITY  RESPONSE 

PLOTS 


1.  SHOT  2  (500  MSEC) 

The  following  vertical  velocity  plots  are  from  the  Shot  2  damping  coefficient 
comparison.  The  Russell’s  error  factor  following  each  of  the  figure  captions  are  only  for 
the  correlation  between  the  simulations  using  the  NPS  damping  values  (a  =  19.2,  and  p  = 
2.09E-06)  and  the  measured  data  from  the  ship  shock  trial  conducted  at  sea. 


DDG-81  SHOT  2 


Grid  1 202  1  7-vz  (V2000V) 


Keel  (x=5328  y— O  z=82) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09 K- Ob) 
NTS  Simulation  (4%  Damping) 


Figure  127.  Keel  Sensor  V2000V:  (RM  =  0.3169,  RP  =  0.2887,  RC  =  0.3800) 
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DDG-81  SHOT  2 

Grid  1 42489-vz  (V2002V) 


Keel  (x=4656y=24  z=85) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09  E- Ob) 
NTS  Simulation  (4%  Damping) 


Figure  128.  Keel  Sensor  V2002V:  (RM  =  0.0679,  RP  =  0.2175,  RC  =  0.2019) 

DDG-81  SHOT  2 

Grid  210430  -vz  (V2007V) 


Keel  (x — 4080  y— O  z=82) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NTS  Simulation  (4%  Damping) 


Figure  129.  Keel  Sensor  V2007V:  (RM  =  0.0879,  RP  =  0.2164,  RC  =  0.2070) 
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DDG-8 1  SHOT  2 

Grid  210808-vz  (V2009VI) 


Bulkhead  (x=4()8()  y=-  1  74  z=177) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NPS  Simulation  (4%  Damping) 


Figure  130.  Bulkhead  Sensor  V2009VI:  (RM  =  0.0882,  RP  =  0.2272,  RC  =  0.2160) 

DDG-8 1  SHOT  2 

Grid  210894-vz  (V2008VI) 


Bulkhead  (x  4080  y=  1  74  z=177) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

NPS  Simulation  (4 °/o  Damping) 


Figure  131.  Bulkhead  Sensor  V2008VI:  (RM  =  0.1200,  RP  =  0.1932,  RC  =  0.2016) 
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DDG-81  SHOT  2 

Grid  220589-vz  (V2010V) 


Keel  (x=3504  y=0  z=86) 


Time  (msec) 


Ship  Shock  Trial  IVIeas tired  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NTS  Simulation  (4%  Damping) 


Figure  132.  Keel  Sensor  V2010V:  (RM  =  0.0827,  RP  =  0.2070,  RC  =  0.1975) 

DDG-81  SHOT  2 

Grid  221  102-vz  (V2012VI) 


Bulkhead  (x=3504  y=-2  1  6  z=177) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

NTS  Simulation  (4 °/o  Damping) 


Figure  133.  Bulkhead  Sensor  V2012VI:  (RM  =  0.1299,  RP  =  0.2211,  RC  =  0.2273) 
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DOG-81  SHOT  2 

Grid  221  188-vz  (V201  1 VI) 


Bulkhead  (x=35Q4  y— 2  1  6  z=177) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NIPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NPS  Simulation  (4%  Damping) 


Figure  134.  Bulkhead  Sensor  V2011VI:  (RM  =  0.0411,  RP  =  0.2240,  RC  =  0.2018) 

DDG-81  SHOT  2 

Grid  2121 96-vz  (V2 1 08 V) 


Bulkhead  (x  4080  y=3  1  2  z=390) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NPS  Simulation  (4%  Damping) 


Figure  135.  Bulkhead  Sensor  V2108V:  (RM  =  0.0809,  RP  =  0.1858,  RC  =  0.1796) 
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DDG-81  SHOT  2 

Grid  221601-vz  (V2013V) 


Bulkhead  0=3504  y  O  z=280) 


Time  (msec) 


Ship  Shook  Trial  Measured  Data 

NJ  PS  Simulation  (alpha  =  19.2,  beta  =  2.09E-6) 

NJ  PS  Simulation  (4 °A>  Damping) 


Figure  136.  Bulkhead  Sensor  V2013V:  (RM  =  0.1049,  RP  =  0.1420,  RC  =  0.1565) 

DOG-81  SHOT  2 

Grid  222060-vz  (V2  1 24V) 


Bulkhead  (x=3504  y=-375  2=390) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NPS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

NPS  Simulation  (4%  Damping) 


Figure  137.  Bulkhead  Sensor  V2124V:  (RM  =  0.1793,  RP  =  0.2311,  RC  =  0.2456) 
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DDG-8 1  SHOT  2 

Grid  222436-vz  (V2125V) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NTS  Simulation  (4%  Damping) 


Figure  138.  Bulkhead  Sensor  V2125V:  (RM  =  0.0214,  RP  =  0.1914,  RC  =  0.1707) 

DDG-8 1  SHOT  2 

Grid  230461-vz  (V2014V) 


Keel  0=2952  y— O  z— 86) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NTS  Simulation  (4%  Damping) 


Figure  139.  Keel  Sensor  V2014V:  (RM  =  0.0590,  RP  =  0.2126,  RC  =  0.1956) 
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DDG-8 1  SHOT  2 

Grid  242399-vz  (V2016V) 


Keel  (x=2544  y=0  z=l  16) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09  E- 06) 
NTS  Simulation  (4%  Damping) 


Figure  140.  Keel  Sensor  V2016V:  (RM  =  0.0169,  RP  =  0.2038,  RC  =  0.1812) 

DDG-8 1  SHOT  2 

Grid  312302-vz  (V2026V) 


Keel  (x=1992  y  O  z=55) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-6) 

NTS  Simulation  (4 °/o  Damping) 


Figure  141.  Keel  Sensor  V2026V:  (RM  =  0.0751,  RP  =  0.2185,  RC  =  0.2047) 
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DDG-8 1  SHOT  2 

Grid  330759-vz  (V2034V) 


Keel  (x=1152  y=-135  z=193) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09  K- Ob) 
NTS  Simulation  (4%  Damping) 


Figure  142.  Keel  Sensor  V2034V:  (RM  =  0.1442,  RP  =  0.1874,  RC  =  0.2095) 
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■  Ship  Shock  Trial  Measured  Data 

■NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

■  NTS  Simulation  (4%  Damping) 


Figure  143.  Keel  Sensor  V2035V:  (RM  =  0.0009,  RP  =  0.1692,  RC  =  0.1500) 
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DDG-8 1  SHOT  2 

Grid  340167-vz  (V2019V) 


Keel  (x=672  y=0  z=197) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09  E- Ob) 
NTS  Simulation  (4%  Damping) 


Figure  144.  Keel  Sensor  V2019V:  (RM  =  0.1327,  RP  =  0.2391,  RC  =  0.2423) 

DDG-8 1  SHOT  2 

Grid  350052-vz  (V2020V) 


Keel  (x=288  y=0  z=21 1) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-6) 

NTS  Simulation  (4 °/o  Damping) 


Figure  145.  Keel  Sensor  V2020V:  (RM  =  0.2477,  RP  =  0.2657,  RC  =  0.3219) 
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2.  SHOT  2  (250  MSEC) 

The  following  vertical  velocity  plots  are  the  supporting  plots  for  the  comparison 
of  damping  coefficients  at  selected  sensor  locations.  These  are  a  subset  of  the  sensors 
listed  in  Section  1  of  this  appendix.  These  sensor  locations  were  specifically  chosen  for 
further  analysis  since  the  simulations  conducted  using  the  NPS  damping  values  provided 
poorer  correlation  with  the  measured  ship  shock  trial  data  than  did  the  simulations 
performed  using  the  fixed  4%  damping  values. 


DDG-81  SHOT  2 

Grid  120217-vz  (V2000V) 


Keel  (x=5328  y=()  z=82) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2. 09  K- Ob) 
NTS  Simulation  (4%  Damping) 


Figure  146.  Keel  Sensor  V2000V:  (RM  =  0.2817,  RP  =  0.2599,  RC  =  0.3397) 
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DDG-8 1  SHOT  2 

Grid  210808-vz  (V2009VI) 


Bulkhead  (  x=4()8()  y=-  1  74  z=177) 


Time  (msec) 

- Ship  Shock  Trial  Measured  Data 

- NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 

- NTS  Simulation  (4%  Damping) 

Figure  147.  Bulkhead  Sensor  V2009VI:  (RM  =  0.0630,  RP  =  0.1868,  RC  =  0.1747) 

SHOT  2 

Grid  221601  -vz  (V20 1 3  V) 


Bulkhead  (x=3504  y=Q  z=28Q) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-6) 

NTS  Simulation  (4%  Damping) 


Figure  148.  Bulkhead  Sensor  V2013V:  (RM  =  0.0700,  RP  =  0.0934,  RC  =  0.1034) 
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DDG-81  SHOT  2 


Grid  222060-vz  (V2  1 24V) 

Bulkhead  (x=35()4  y=-375  z=39Q) 


Time  (msec) 


Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NTS  Simulation  (4%  Damping) 


Figure  149.  Bulkhead  Sensor  V2124V:  (RM  =  0.1418,  RP  =  0.1720,  RC  =  0.1976) 

DDG-81  SHOT  2 

Grid  330759-vz  (V2034V) 


Keel  (x=l  152  y=-135  z=  I  93) 


Time  (msec) 

Ship  Shock  Trial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09E-06) 
NTS  Simulation  (4%  Damping) 


Figure  150.  Keel  Sensor  V2034V:  (RM  =  0.1138,  RP  =  0.1456,  RC  =  0.1638) 
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Vertical  Velocity  (ft/s) 


DDG-8 1  SHOT  2 

Grid  350052-vz  (V2020V) 


Keel  (x=288  y=0  z=21 1) 


Time  (msec) 

Ship  Shock  X rial  Measured  Data 

NTS  Simulation  (alpha  =  19.2,  beta  =  2.09H-6) 

NPS  Simulation  (4%  Damping) 


Figure  151.  Keel  Sensor  V2020V:  (RM  =  0.2161,  RP  =  0.2423,  RC  =  0.2877) 
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